{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/imarkelo/anaconda3/lib/python3.6/site-packages/cmocean/tools.py:76: MatplotlibDeprecationWarning: The is_string_like function was deprecated in version 2.1.\n",
      "  if not mpl.cbook.is_string_like(rgbin[0]):\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "'2019-01-12 12:45:42.355345'"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "\n",
    "from matplotlib import rc\n",
    "import matplotlib.ticker as tkr\n",
    "import matplotlib.dates as mdates\n",
    "from datetime import datetime, timedelta\n",
    "from matplotlib.colors import ListedColormap\n",
    "import seaborn as sns\n",
    "sns.set_style(\"whitegrid\")\n",
    "sns.set_style(\"ticks\")\n",
    "import scipy.io as sio\n",
    "import cmocean\n",
    "\n",
    "import sys\n",
    "sys.path.append('../')\n",
    "\n",
    "\n",
    "import warnings\n",
    "warnings.filterwarnings('ignore')\n",
    "\n",
    "%matplotlib inline\n",
    "\n",
    "\n",
    "sns.set_style(\"whitegrid\")\n",
    "sns.set_style(\"ticks\")\n",
    "\n",
    "rc('text', usetex=False)\n",
    "rc(\"savefig\", dpi=90)\n",
    "rc(\"figure\", dpi=90)\n",
    "\n",
    "plt.rcParams['figure.figsize'] = 6, 4\n",
    "\n",
    "import pandas as pd\n",
    "\n",
    "pd.options.display.max_columns = 999\n",
    "pd.options.display.max_rows = 400\n",
    "\n",
    "import h5py\n",
    "\n",
    "from datetime import datetime\n",
    "str(datetime.now())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import sys\n",
    "# sys.path.append('../')\n",
    "from porousmedialab.column import Column\n",
    "from scipy import special"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Simulation started:\n",
      "\t 2019-01-12 12:47:27\n",
      "\n",
      "\n",
      "Estimated time of the code execution:\n",
      "\t 0h:00m:02s\n",
      "Will finish approx.:\n",
      "\t 2019-01-12 12:47:30\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhMAAAFfCAYAAAD00U6uAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAN1wAADdcBQiibeAAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XdYlfX/x/HnYQ9RxBUKlisRt7hzfNPUDA1HZmpLrczV+LWob7ZL0oZly10aWpm4Tc2tuQVXaimpiOJAhsg4rPP7wy8nUNSDwDkIr8d1eV2c+9z3fd58tHhx35/78zaYTCYTIiIiIrfIztYFiIiIyO1NYUJEREQKRWFCRERECkVhQkRERApFYUJEREQKRWFCRERECkVhQkRERApFYUJEREQKRWFCRERECkVhQkRERApFYUJEREQKRWFCRERECsXB1gUUtT179ti6BBERkdtCQEBAkZyn1IUJKLrBkRszGo0AODs727iSskHjbV0ab+vSeFuX0Wjk4MGDRXY+3eYQERGRQlGYEBERkUJRmBAREZFCKZVzJkRESguTyURmZiYmk8nWpRSrjIwMAAwGg40rKR0MBgMODg5WG0+FCRGREio1NZVTp06RlZVl61KKXU5YUpgoOvb29vj6+uLq6lrsn6UwISJSAplMJk6dOkW5cuWoXLmyrcspdgoTRS82NpZTp05Rr169Yh9XhQkRkRIoMzOTrKwsKleujJOTk63LKXbZ2dkA2NlpKl9RqVy5MomJiWRmZuLo6Fisn6W/NRGREqi0z5EQ67HGvyWFCRERuS2ZTCbOnDlj6zIKJSUlhbi4OFuXUWg2CROzZs0iODg43/fS0tJ46aWXaNmyJZ07d2bhwoVWrk5ERG7FI488QpcuXcy3LAqjS5cu7Nix44b7fPzxxyxYsACA3bt3ExgYeMuft2PHDrp06XLLx1/tscceIyws7Kb7Pfrooxw9ehSA7777jrfeeqvIarAmq4aJ9PR0Jk2axIQJE667z2effUZqaiqbN2/mq6++IiQkhOPHjxfoc07FpTBizm6embObU3EpnIpLYfTccGISUhk9N9y8Lfc+t3pcfvsAJKVl5HktIlKa/fPPP8TFxVG1alU2btxolc9MSEgwf92yZUuWL19ulc8tSrm/h2effZb33nvPhtXcOqtOwAwODiY5OZkBAwaQnp6e7z7Lly/nu+++w9XVlcaNG9OrVy8WLFjAyy+/bPHnfLD8EKv+PGd+7Whvx/L9Maw/cp6U9CuPWGVmZefZZ+pjLW/puPz2GXtvXR6Ztp2ElCvPTX89uAWn4lL4YPkhTMC4QH98vdyAK6EjOOwAwff7mbeJiNxuFixYwL333stdd93FTz/9xL333gtc+Q29VatWrFixgtjYWO69914+/PBDnJycOHHiBB988AFHjhwhOTmZtm3bMnHiRMqVK2c+75dffsnRo0eZPHkycOVx2XvuuYdhw4axdOlSDAYD8fHx9OjRg9dff51169YB8P333zNr1iySk5Np3749H330EeXKlWPhwoXMmjWL6Oho3NzceOaZZ3j88cdv+L2FhoYybdo0UlNT8ff357333sPX15dLly4xfvx41q9fj5OTE0FBQTz33HPXTHasX78+a9euxcfHB7hy1WX8+PH8+uuvnDlzhqeffpqQkBAiIyM5ffo0ISEhNzx3cHAwFSpUYNeuXZw4cYJWrVrx8ccf4+npWWR/nwVl1SsTwcHBTJkyhSpVquT7fmJiIrGxsdSuXdu8rVatWhw7dqxAnxMdewlTViam7Cxi4lN4oWtd3JzsSUnPws3Jnhe61iMmMc28/9nENGISU6/Z9uJ99W56XH775AQJNyd73nygAfBvwFn95zneX37IfIVjyPQdLN8fQ8jKI4CuaIjI7ScrK4vFixfTt29fevXqxc6dOzl9+rT5/dWrVzN79myWLVvGzp07WblyJQDjxo0jICCATZs2sXr1ak6ePMmyZcvynDswMJDNmzeTmpoKwIYNG2jYsCFjxoyhd+/ejBgx4ppbA2vXrmX69OlMmzaNP/74g+zsbCZNmkRUVBQffvghn376KeHh4bz//vtMmDCBy5cvX/d7u3jxIp9//jnz5s1j69at+Pr6MmXKFHP9CQkJ/P777/z6669s27bN/J4lJk6cSPXq1Zk2bRoPPPBAnvdudu5ly5bx2WefsX79es6fP8+8efMs/tziYNUrE1WrVr3h+zn/WHIvsOHi4kJaWtr1DsnXiv/rav46Clj2/P9eGK5kp7s/tAODATs7AwaDHVEmqPGSATDgYG8HBjtOmWD5a4Yr+zk4gb0jDb9zxODgSDlXV+wcnFiTZaDpZEcMjs44uXmQ4OhOmz8WYOfsjmu58sQ5uPHop3F8New/nIlPNdcUE5/KW4sOsP7vWAAaVy/Pe4H1OXAylse+DychNYPsrGwmPdwYgMtpmby55DAvd6uLT8XiX3zEUte7uiTFQ+NtXbYe74yMDEwmEyaTyTwHYejQoSxevLjYPzsoKIhZs2ZZvP+GDRu44447uPvuuwHo2rUrP//8My+88AImk4mgoCDzWhmtWrXi5MmTZGdnM378eCpXrkxKSgpnz56lQoUKXLhwgezsbPP3XqtWLWrWrMnGjRvp3r07K1eupGfPnnn2ufrr3377jYceeoi6desC8NZbb3H58mWqVq3KokWLqF69OufOncPOzo6MjAzi4+OvGescTk5OZGRk8PPPP9OzZ0/eeust7OzsSElJYc2aNSxZsgR3d3fc3d157rnneP/99xk1atQ157v66/zqzvnaknN369aNmjVrAtCxY0eioqKuqT3nfOnp6dc80VHU/75L1DoTLi4uwJVJmG5ubtd8bSmHitXx8XQhOzubMwmpVwbRZKKSuyMXL18JJiZTNhVdHDBhIiE5Hbgy6C4OdmRnZ5OSkQEmMJmyyM7KzHP+/ObdJl+nlnOA/wdgcHanXMWquHhWYaNbFRwr+eJQyRfHSr6YTB5cNmaag4Sbkz3BPeoBV4LE8DkR7I1OBDAHDBGRkiIsLIxjx47RsWNH4Movhq6urowePRoALy8v87729vbmH2zHjh1jxIgRJCYm4ufnx+XLl/N9jLFnz56sXr2ajh07smXLFsaNG3fDeuLi4ggICDC/rly5MpUrVyY7O5vQ0FAWLVpExYoVadq06U2/N3d3d6ZOncrUqVOZOXMm1apV480336R+/fpkZmZSo0YN877Vq1fn7NmzNz3nzVy6dOmm5849pg4ODkUy6bUwSlSY8PT0xMvLixMnTuDv7w/A8ePHqVWrVoHOU+OZqTT2qQCA4X8/hAHcnOxxTf93Wdom/9tnf6598tvm6mhHSqoRU1Y6pswM6ldxITvTyJHT8Zgy0zFlpOGYmUpqchLZaclkp12mklMGjSrbs3rXETKTLpKVFEvS2eMknb12MmmMmycNvqqPc/X6lL/Tn+Sqd/PCrwcZF9iAYT/sJiElg+Y1PZkwoCnOzo4lbp6Fs7OzrUsoUzTe1mWr8TYYDOY/OQs5/fDDDzap5Ubi4uLYvHkzYWFhVKhQwbz94YcfZu3atdd8DzmvMzMzefHFF5k8eTIdOnQgOzubMWPGmPfNfVyvXr3o27cvmzdvplmzZuarHLn3yf11tWrVuHDhgvkzjx07xpYtW6hUqRJbtmzht99+w9PTk0uXLhEWFpbvWOdISEjAycmJGTNmkJqaSmhoKK+++irbtm3D0dGRmJgY88+o06dPU6lSpWvqMRgMZGdnY2dnh8lk4tKlS/nWnfN11apVLT731eOQW852JyenYl/4rESFCbhyf+zLL7/kk08+4cSJEyxbtozvv/++QOdo4lPBHAaa+FTgs4eb8uBXf5jnNSwefQ8vzd+XZ58pjwUwYs6e6x7n7ubC4tFdcx3nQsvmNfPsU/V/5351wX4iohI46mRPtTpXwovJZMK/sgOp8ef58/ARMi6eIuPiKRwSz3Ap5h9Sj+0g9dgOEgB7R2cu+DZm48IA3Oq1pXzlO/hmcAs8XK4Eicdn7iQi6soM4K8HtyiScRcRuRWLFy+mefPm5lsKOXr16nXD+/jp6ekYjUbc3NwwmUysX7+ezZs3U69evWv29fX15c477+Trr7/miSeeMG93cnLKd75DYGAgb7zxBr169aJ69ep8+eWXVK9eHUdHR/Of5ORkJk6cCFxZbfR64uLiGD58OHPmzKFhw4aUL1+eChUqYGdnR2BgIBMmTGDChAmkpqYyefLka+Y+ANx5550sW7aMESNG8P3335Oc/O+1bEdHx2u+h4Kcu6QoEYtWNW/enN27dwPwf//3f1SsWJGuXbsyZswYXn/9dRo0aFCg83lXcDF/fUcFFz5fc9QcJFLSs5i09ug1+3hXcL2l4/Lb55vBLcyvPd0c2fZ6F5r6enL4YhYnsivRrusDHF46la4jP6Dio5/j+8IvVBsygYr3Dse7aUdcnJ1I/Wc38WumcPrbYfzz/Ss8MOptDp88w8Cp24iISqCJTwVC+v17y0MTN0XEFsLCwujZs+c124OCgti5cyf//PNPvseVK1eO//73v4wZM4Y2bdowa9Ys+vbtS2RkZL77BwYGcvz4cbp3727e1r17dxYsWMCrr76aZ98OHTowYsQInn76aTp16oSLiwvPP/88ffv2pVq1anTs2JEePXrg4uJC3bp1r/uZALVr1yY4OJixY8fSvHlzfvrpJz755BMA3nzzTTw9PenRowcPPvggAQEBvPDCC9ec45133mHZsmW0adOGqKgoWrT495fAoKAgXnrpJX7++ec8x1h67pLCYCpla7bu2bOHqrUa5HkMEyBk5RHefKABH6w4TPD9fgDXPKp59eOblhyX3z4hK4+wfH8Mnm6O/PR0W/y8yzNizm7zI6XdG1bj0wFNGTJ9h/lKSBe/qsSnpBMRlYCrvYn4E3+SemwnyYc3k5V0AQCDozPuDbtQPuBBApo1IvSpNtdcrQhs4m21qxVGoxHQZXdr0Xhbl63HOz09ncjISOrUqaPeHMCKFStYsmQJ3333nTXLuq3d6N+Q0Wjk4MGDeeaWFEapDBNFNTi36lRcCiErjxDSrzEeLo7mbblDSE7gaOJTgTsquDAu0B8HOwNdP9tovqLx2/MdGTlnN9u3bSX5z/UkH9qAKePK/+Bc67Wl7YBnCX31YUaFhhMRlUDzmp7MHtba/JnFzdb/sy1rNN7WZevxVpi4IiUlhaioKN555x2GDh1Kjx49bFHebcmaYaLEzZkoDXy93K65OuDr5caUx1qaX+dc5cgdOEbPDTcHiZ+ebot3BVe+fbQlXc8n4+LbCM//DOXyvtWk7V1G6tHtrP9oB3UXzqRi5yep7F3DPK8iR0mbqCkiUlBxcXEMGjSIjh075rnFISWLwoSN5Bc4rg4YSWkZjMoVMGaM7MYHy30ID+hNuRMbOLLyB1IObyL16A6S73mEZz2d+PGZezRRU0RKDR8fHyIiImxdhtxEiZiAKVfkBIycqwvBYQfMty82v3ovAXd6MXtYa5reVZnkut2oMWI6Fdo/gsmURcLGH/jt/Sd48N1QYhJTzUGieU3PPBM1RUREiprCRAkWfL8fgU28rzsPovGdVTm8YiZd/zsblzubkBEbxYaJT9Ow33OEn4yz+hwKEREpmxQmSrCrr1TAlasV+6MTaV7Tk5+euTKvYuHrD9Hj5a+o2PVpAOLXTefCgvd5q/ud18yh0OOjIiJS1DRn4jaT38RNDxdHvn20JT3jUnGp2ZgLiyeQGrmLtu3uYePvv9G2mb/mUIiISLHRlYnbTH5XK3ImaiakZNCmZQt27tiGe92WpMdG0aF9O35a+rvmUIiISLFRmCgFck/UnD2sNc3r+nBk+3q82vQhK/USg/v1ZuvG9ZpDISJFrn79+uYVIXNER0dTv359G1V0pUtoYRa3Cg4OZvLkyUVYUemnMFEK5DdR06dSObYvnk3Fe4djyjRyfsG7XP57+zXHah6FiBTWrFmz2L9/v63LMHvvvfd49tlnbV1GmaIwUQpc79bHqwv2U751X7x6jIasLNZMDqbbS1+SlJZh3ufxmTtZvj+GkJVHbFW+iBShU3EpjJ4bbv7vPEdx/uLw4IMP8sYbb5Cenn7Nezt27KBLly7m17mvWuzYsYOBAwfyzjvv0LJlS+6//352797N6NGjad68OQMHDuTChSvtBJKTkxk3bhz33HMP9957L9OnTzefs0uXLvz3v/+ldevWTJ8+Pc+Vhfj4eJ5//nkCAgLo3Lkzc+fOBa6srPnGG2/QpUsXmjRpQr9+/fjrr7+KfGzKCoWJUij3ZMvmNT05+PNEmg5+FbIz2Tn1v9z/+jStRSFSSuUs1f/4zJ1W+8VhxIgR2Nvb8/XXXxf42L179+Ln58fOnTtp3rw5Tz75JIMHD2bbtm04OTmZO4+GhIQQFxfHqlWrmDNnDgsWLGDFihXm8yQlJbFlyxYeeeSRPOcfN24cDg4ObN68mVmzZvHFF1/w559/MmPGDOLi4li+fDm7du2iVq1aTJkypXADUYbpaY5S6Oo5FB4ujmye8T5tUy9zaOHXbPv2VQIumXCqcpfmUYiUMiH9GnMmIZWIqAQen7mTb4a0yNO/pzh+cXBwcGD8+PEMGjSI7t27U6FCBYuPdXd3Z+DAgZhMJgICAvjrr7+45557AGjZsiUxMTGYTCYWL17MkiVLKFeuHOXKlePJJ58kLCzM3Ja7e/fuODk55elBYTQa2bBhA7/99htubm7Url2b2bNn4+3tzaOPPorBYMDBwYHo6GjKly/PsWPHinZgyhCFiVLoeo+Pbp/7OR0GpLB/2SzOL3gf78c+5ZshXRQkREoRDxdHZg9rbb7y2G78OoBi/8XB39+foUOH8vrrr/Pll19afFyFChUwGAyYTCbs7e3x8PAwv2dnZ4fJZCIuLg6j0chDDz1kfi87OxsfHx/z68qVK19z7sTERDIyMqhWrZp5W84tlpMnT/L2229z+PBhateujaura4G+X8lLtzlKofzmUOSodf8w3Pw6kpV4jgsLP2LE99uvubcqIrc3DxdHvhmSdy2Zb4bk//+EojRq1Ciys7PzPElhZ2dHZmam+XVCQkKeYwwGw03P6+npiaOjI6tWrWL37t3s3r2btWvXMm3atBuep1KlSjg6OnLu3DnztgULFhAeHs67775Ls2bN2LZtG/PmzaNTp04F+l4lL4WJMiLnnuneU4l0HfEODZu2wHj6EOtnfXzNvVU93SFye0tKy2BUaHiebaNCr52UWdScnJz46KOPWLJkiXmbr68vsbGxbNu2DaPRyNSpUwt8Xnt7e3r27Mknn3xCSkoKiYmJPPfcczd9/NPe3p4ePXowefJk0tLSiIyM5NNPP8XV1ZVLly7h6uqKnZ0dx44dIzQ0lIwM/WJ1qxQmyojc8yhCn+1Iy6c/ws7dk8v7VrF5RRjBYQf0dIdIKXD1BOxtr3eheU1P8xyK4g4UTZo0YdiwYebXd9xxB6+88gqvvPIKXbt2pV27drd03rfeeguAbt260b17d2rUqMGrr7560+PGjRtHVlYWnTt35umnn+bVV1+lQYMGBAcHs3DhQpo3b87zzz9PUFAQJ0+ezHMVRSxnMJlMJlsXUZT27NlDQECArcsocU7FpRCy8kie9uY9g7/jjy9fwM7BkZ9WrGPu36ZrJm7eiNFoBMDZ2dka30KZp/G2LluPd3p6OpGRkdSpUyfPpMKbGT03nOX7Y/L8d5w7YAQ28S6Ry+lnZ2cDV26LSNG40b8ho9HIwYMHi+znpf7Wyoir51F4uDjyW8izNOj1FNkZRoYMGkT4P+f1dIfIbS6/RexyJmUGNvE2T9AWKUoKE2WYh4sjv//wGS53NiMjNoqEzT9aZZKWiBSf603A9nBx5OvBLfD1crNRZVKaKUyUYUlpGYyZt5dKDzyPwdmdSzsXMvDdWXq6Q0RECkRhoozKfQ+1VaN6TJo0CTCxbdb7DP5mgwKFiIhYTGGijLp6lcyxI4bTtUcgmQln2Th3MsFhB8z76nFRERG5EYWJMurqSVoGgwHfoOexc3YnKXw5XSpdBtQMTMRWLFnMScQS1vi3pOW0y6icSVq5fTn0P0RsGsm+nz7hqRHP0iZ8J2N/2qdmYCI24ODggL29PbGxsfkuFV3a5KxSoBBVdGJjY7G3t8fBofh/1CtMiJmHiyMbpr9P7W3LiT95mKaDXsOj+QN6XFTEBgwGA76+vpw6dYrExERbl1PsFCaKnr29Pb6+vroyIdbn6e7CvB+mc/+9HUnY+ANud7dXMzARG3F1daVevXpkZmZSytYXvEZ6ejpAgRbokuvL6YhqrXCmMCF5JKVl8N2fJjxaBJK0ZykJW0IZFXqXrkyI2IjBYMDRsfT/t5cTlhQmbk+agClmuR8X7ThwJBUqeJK8fzU7wvdbZU1/ERG5PSlMiFnux0V/Gnsfb701DlN2Ftk7fiQiKiHP46IiIiI5FCbE7OrHRUePHk2tWrU4d3ArTe1PaU1/ERHJl8KEmF29pr+zszMvv/kuACdXTKGGp4t536S0DF745QDR8ak2qVVEREoOhQm5oYPO/jh512ff3r2E/vQL8O/cit/+PMcnvx+zcYUiImJrChNyQx/3b0LLh54FYPTLb3A6Ptk8SbOZTwU+eLCBjSsUERFbU5iQG/JwceS3j0dRqU4TkmKO0+zJ98yTNGc81pxyLnq6WESkrFOYkJsq7+rEd59/DEDClrmYsrP4ZkgLBQkREQEUJsQCSWkZhJ7ywNm3EZlx0SQf3sSo0HAup2XaujQRESkBFCbkhnIvZNWq/wgAjLvmE34yjuFzIhQoREREYUJuLPdCVivGj6BTp05cPhfFHfEH2BudyJtLDtu6RBERsTGFCbmhqxeyeu211wDI3ruY+/2r8nK3ujauUEREbE1hQm7o6oWsevbsSePGjdm9ayf9vRPwqehq4wpFRMTWFCakQKLjU6na4WEAPv30U/P2pLQMRs8N51Rciq1KExERG7F6mNi3bx99+vShWbNmDB48mKioqGv2ycjIYNy4cbRt25Z27doxfvx4srOzrV2q5CNk5RGOujfGtWI1Vq1axYEDB8yTNJfvjyFk5RFblygiIlZm1TBhNBoZPXo0w4cPZ+fOnbRv357g4OBr9gsNDeXMmTOsW7eO5cuXs2XLFhYtWmTNUuU6Qvo1pkWtyri0eBCA8RM/NT/t0bymJyH9Gtu4QhERsTarhont27fj6elJ7969cXJyYuTIkRw9epTIyMg8+0VFRZGVlWW+GmFnZ4eLi0t+pxQr83BxZPaw1tzzwAAMTm6ELfiV3YdP0Lymp3mSpoiIlC1WDRPHjx+ndu3a5tf29vb4+vpeEyYGDBjAX3/9RatWrWjXrh21a9fmgQcesGapcgMeLo5MHd6Bck26QVYmSftW8s2QFgoSIiJllFXXQ05JSbnmCoOrqyupqXnbWKenp9O7d29Gjx5NUlISzz77LKGhoQwZMsSizzEajUVWs1zrclomI+dE4NEikKTdS7gcsYIR3+9g1pOttMR2MUpPT7d1CWWKxtu6NN7WVdTjbdUrE66urqSlpeXZlpqairu7e55tr7/+Or169aJChQr4+Pjw7LPP8uuvv1qzVLmOy2mZDJ8Twd7oRAIa1ec/Xe8j63Ic29b9phUxRUTKKKv+Glm7dm0WLlxofp2VlUVUVBS1atXKs9/Zs2fJyMgwv3ZwcMDBwfJSnZ2dC1+s5Ov/FvzJ3uhEmtf0ZPqQZmyr9Rwb1v5O9sHf2NugE28t/4uvB7ewdZmlmv59W5fG27o03rcnq16ZaNOmDRcvXmTRokWkp6fz7bffUrNmTerUqZNnv44dO/LFF1+QlJTEhQsXmDZtGj169LBmqXIduVfELOfiwH333Yefnx9x/xyglUcCwff72bpEERGxMquGCRcXF6ZMmcKcOXNo06YNW7duZdKkSQAEBgayZMkSAN59912qVq1Kt27d6Nu3L+3bt+fJJ5+0ZqlyHVeviGkwGBgzZgwA7kd/x9fLzZbliYiIDRhMJpPJ1kUUpT179hAQEGDrMsqEnImu6enp1KhRg4yMDM6cOUPFihVtXFnplDPeugxsHRpv69J4W5fRaOTgwYNF9vNSy2lLoSVk2OPbqjtpaWn8+OOP5u1aYltEpGxQmJBCC1l5hDjfjgBMmToVk8mkJbZFRMoQLQoghRbSrzFnElJZsbIOfx48yPK1m5jxt4OW2BYRKSN0ZUIKLWeJ7Qb39gXgkf/70BwktMS2iEjppzAhRcLDxZFfJryMwdGZlCObyDamaIltEZEyQmFCikRSWgavLT2Gm18nTBlGkg9tYFRoOElpGTc/WEREbmsKE1JoOZMtI6ISCOjRH4Csw2uIiErg8Zk7FShEREo5hQkptOCwA+Y5EkvfG0qjRo1IPPU3te0uEBGVQHDYAVuXKCIixUhhQgot9xLb5V2dzKuVNkzeR2ATby2xLSJSymkFTLll11uxLiYmBh8fH6pUqUJ0dHSBmrTJ9WmFQOvSeFuXxtu6tAKmlHje3t706NGDc+fOsXr1aluXIyIixUxhQorFE088AcAPP/xg40pERKS4KUxIkTsVl8K6VF8qVKjA4sWLiY+PB9SrQ0SktFKYkCIXsvIIq47E49moE0ajkfnz56tXh4hIKaYwIUUupF9jmtf0JO2uK82/ps2cZV6HQr06RERKH4UJKXI5vTratmuHQ0Vvdu/Yzs59h9SrQ0SklFKYkGLh4eLIt48G4O5/LwDJhzaqV4eISCmlMCHFIiktg1Gh4bj7dwYg5fAmRv64R0tri4iUQgoTUuRy9+po3awhjZs2I+PiKXbuiVCvDhGRUkhhQopc7l4ds4e15vFHhwBQLnqHenWIiJRCChNS5HL36vBwcWTgwIEAGP/awgONqqlXh4hIKaMwIUXO18uNrwf/O9nS19eXjh07EhV1kiG1jPh6udm4QhERKUoKE2IVgwYNAmDevHk2rkRERIqawoRYxUMPPYS9vT3z588nMzPT1uWIiEgRUpgQq6hSpQrdu3fn/PnzrFu3ztbliIhIEVKYEKs4FZdC+p3tgLy3OtT8S0Tk9qcwIVYRsvIIf7v4YefoxMKFCzEajWqUgRqOAAAgAElEQVT+JSJSSlgUJiIjI/n999/NrzMzM3njjTc4duxYsRUmpUtIv8YE1KuOc60AEhMTmb9khZp/iYiUEjcNE/v27WPgwIH89ddf5m2pqalkZGTwyCOPcPDgwWItUEqHnOZf/u27AfDse9/kWdhKPTtERG5fNw0Tn3/+Oa+99hpjxowxb/Pw8GDixImMGTOGzz//vFgLlNLDw8WRue+NAntHUo9uw5SVoeZfIiKlwE3DxOHDh+nXr1++7w0aNEhXJsRiSWkZvLb0GK61WpBtTCbt5H5GhYarV4eIyG3upmHCzs7uuusCGAyGIi9ISqfczb8a/O9Wh+vpXUREJaj5l4jIbe6mYSIgIICFCxfm+96iRYto0KBBkRclpU/u5l+LPn4OR0dHLh3ZStMa5dT8S0TkNudwsx3GjBnDkCFDiI6O5r777sPLy4uLFy+ydu1a5s2bx9SpU61Rp9zmcpp7hfRrjIeLI926dWPFihUMq53C6kreav4lInIbu2mY8PPz4/vvv2fChAnMmDEDk8mEg4MDrVu3ZtasWTRp0sQadcptLqf5V44BAwawYsUKVixZxJQpU2xYmYiIFJbBZDKZLN05PT2dxMREPD09cXQsmTPw9+zZQ0BAgK3LKBOMRiMAzs7OBT42Li6OatWq4enpSUxMDA4ON821ZV5hxlsKTuNtXRpv6zIajRw8eLDIfl4WaAVMJycnqlSpUmKDhNw+vLy86Nq1K7GxsWzatMnW5YiISCFoOW2xmQEDBgAwf/58G1ciIiKFoTAhNhMUFIS9vT1hYWFkZWXZuhwREblFChNiE6fiUnh7dRQdO3Xm/PnzbNu2DVAXURGR21Ghw8Q777zDe++9x5Ej6voolgtZeYTl+2OIq9IUgIULF6qLqIjIbarQYcJkMjFmzBjCw8OLoh4pI0L6NaZ5TU8uVroSJhaELeSxGTvURVRE5DZU6DDx7rvv4uXlxeDBg2+67759++jTpw/NmjVj8ODBREVF5bvfjBkz6NChA61ateKll14iLS2tsGVKCZPTRbRVo7o4edfj5Inj7NyzV11ERURuQwUKExERESxevJhFixbl+WMJo9HI6NGjGT58ODt37qR9+/YEBwdfs9+KFSuYN28e8+bNY8OGDcTHxzNt2rSClCm3CQ8XR74Z0gK3eu0ASDm6XV1ERURuQxavFDR+/Hjmz5+Pn58f9vb25u0Gg4E+ffrc9Pjt27fj6elJ7969ARg5ciQ//PADkZGR1KlTx7zf/PnzGTNmDL6+vgBMmDBBVyZKqaS0DEaFhuNWrx0Jm2aT8vc2RoWG68qEiMhtxuIwsXjxYn7++Wfq1at3Sx90/PhxateubX5tb2+Pr6/vNWHi8OHD3HfffQQFBREbG0uPHj3yvYJxIzkrqUnxSk9Pv+VjL6dlMnxOBHujE2nZtCEH69Tln8hj7DrwF49NNzHjseaUc9GqmLkVZryl4DTe1qXxtq6iHm+Lb3O4ubnh7e19yx+UkpKCi4tLnm2urq6kpqbm2Xbp0iXCwsKYMmUKS5cu5dChQ+rdUAq9ueQwe6MTaeZTgZmPt6Bf3ytXtyqe38ve6ETeXHLYxhWKiIilLP7Vb9CgQYwdO5ZHHnkELy+vPO+1atXqpse7urpec7siNTUVd3f3PNscHR157LHHuOOOOwAYNmwY3333HWPHjrW0VK3tbmW3Mt5vBPpjZ29n7iL60EMP8cknn+B5YR+d+j9J8P1++nu8Do2LdWm8rUvjfXuyOEzMmzcPgI8//jjPdoPBwNq1a296fO3atVm4cKH5dVZWFlFRUdSqVSvPfnfddReXL1/Os18BepHJbeLqLqKtWrXC29ubrX9sIWyBL1W83GxYnYiIFITFYWLdunWF+qA2bdpw8eJFFi1axAMPPMDUqVOpWbNmnvkScGWJ5VmzZtG1a1ecnZ2ZOXMm3bt3L9RnS8lnZ2dHnz59+Pbbb1m2bBlDhw61dUkiImKhAj0aunv3bsaNG8fTTz9NcHBwgbo9uri4MGXKFObMmUObNm3YunUrkyZNAiAwMJAlS5YA8OSTT9K/f38GDx5Mjx49aNSoEcOHDy9ImXKbynkqKPcVLBERKfkMJgvvISxZsoQPP/yQgQMHUr16daKjo5k/fz6vvPIKDz30UHHXabE9e/YUWX92ubGcp2aK6h5neno6VatWJS0tjdjYWMqVK1ck5y0tinq85cY03tal8bYuo9HIwYMHi+znpcW3Ob777jumTZtGkyZNzNu6d+/OSy+9VKLChNy+nJycCAwMZO7cuaxatYr+/fvbuiQREbGAxbc5Lly4gL+/f55t/v7+xMfHF3lRUnb17dsXwOKVVUVExPYsDhONGzdm5syZebbNmDGDRo0aFXlRUjadikth9aU7cHZ2ZtmyZWRkZABqSy4iUtJZfJvjv//9L0899RRz587ljjvu4MyZM7i6uvLdd98VZ31ShoSsPMLqvxPxrBfAuYNb2bhxI206dObxmTuJiEoAyPM4qYiIlAwWh4k6deqwatUq9uzZQ1xcHNWqVaNp06Y4OqqHghSNkH6NOZOQyiaflnBwKz/+PJ+v/3ZVW3IRkRLupmFi2bJl9OrV65p72NHR0URHRwNY1OhL5GZy2pIPTIrnt5WT+fGnBdTw6k2LOyuq+ZeISAl20zCxePFievXqxYIFC/J939KuoSKW8HBxZMazXan1VQOM0YdIP3uMb94YoSAhIlKC3TRMTJs2DYDPPvuMKlWqXPN+ZGRk0VclZVZOW3LXum0xRh8i9eh2RoW20pUJEZESzOKnOXr06HHNtuTkZAYMGFCkBUnZlZSWYZ5s2bJzNwBMJ3cREZXA4zN3kpSWYeMKRUQkPze8MnHmzBkCAwNJS0vDZDLRoEGDa/bp0KFDsRUnZUtw2AHzZMvZw7rTdr4/hw4doo1rMhFRV97X0xwiIiXPDcNE9erVWbNmDampqTz22GP8+OOPmEwmDAYDcGXFwvxufYjciuD7/QDMbcn79OnDoUOHuNflJPXqdTW/LyIiJYvFvTnyk5WVRWRkJHfffXdR1lQo6s1hPcW9lv6uXbto3bo1HTt2LFBTudJKvQusS+NtXRpv67JZb46VK1fywQcfcPHixTzby5cvz44dO4qkGJHcAgICqF69On/88QcXLlzQVTARkRLK4gmYEydOZOTIkbz//vs88MADLFy4kM6dO/Pkk08WY3lSltnZ2REUFER2djbLli2zdTkiInIdFoeJuLg4hgwZQseOHTl+/Dh+fn589NFHzJ8/vzjrkzIuKCgIuLLeiYiIlEwWhwlvb28uXLhAtWrVOH36NOnp6Xh4eJCQkFCc9UkZd++991K+fHlWr15NSooafYmIlEQWh4m+ffsycOBALly4QJcuXRg5ciTPPfccDRs2LM76pIxzcnKiZ8+epKamsnr1aluXIyIi+bA4TDz99NOEhIRQvnx53nnnHVq3bk3t2rX57LPPirM+KeNOxaUQV6UpkPdWh9qSi4iUHBY/zdGvXz/mzJljfmxnxIgRxVaUSI6QlUc44lAHg70DS5YuJTMzk9RMk9qSi4iUIBaHiUuXLpGYmIi7u3tx1iOSR05b8gu+jYk7EcHS1euZfcJNbclFREoQi8PE3XffTZ8+fWjSpMk1z/uPHz++yAsTgX/bknfc1IV9JyJ44u2v8er69P+W3FbzLxGRksDiMOHv74+/v39x1iKSLw8XR2a+/SwBP39K6rEdmLo8xTdDWihIiIiUEBaHCVdXV4YPH37N9i+++KJICxK5WlJaBu+vP4eTdz3SY46SEXuSUaHhujIhIlJC3DBMnD9/nq1btwLw5ZdfUqlSpTzvX758mdmzZ/P8888XX4VSpuVuS16n5X84vPQoFc7vJSLqLh6fuVOBQkSkBLhhmPDy8mLdunXEx8eTkZHBggUL8rzv5OTEm2++WawFStmWuy356w+Pos3SabieCcev55NERCWoLbmISAlwwzDh4ODAl19+CcDnn3/Oiy++aJWiRHLkbkteztmBunXrEhEeTuhP3ny1y1VtyUVESgCL50y8+OKLREREEBUVxdVdy/v06VPkhYkA+Hq55bnyEBQUxKeffsq61b/x9ejRNqxMRERyWBwmQkJC+OWXX/Dz88Pe3t683WAwKEyI1fTp04dPP/2URYsWMVphQkSkRLA4TCxatIiff/6ZevXqFWc9IjfUrl07qlSpwoYNG0hISMDT09PWJYmIlHkW9+Zwc3PD29u7OGsRuSl7e3t69+5NZmYmK1assHU5IiJCAcLEoEGDGDt2LKtWrWLXrl15/ohYU1BQEHDlapmIiNiexbc55s2bB8DHH3+cZ7vBYGDt2rVFW5XIDXTr1g03Nzd+++03jEajufmciIjYhsVhYt26dcVZh4hFTsWlELLyCF3v68bSJYtZt24dPXv2JCktg+CwAwTf74evl5utyxQRKVMsvs0BsH//ft59911Gjx5NfHw8M2fOvOYxUZHiFLLyCMv3x3C6QiMAFi9ebF4lc/n+GEJWHrFxhSIiZY/FYWLx4sWMHj2a8uXLs337drKysggLC2PixInFWZ9IHiH9GtO8pifnKzYEgx0LFy3isenb1ZJcRMSGLA4T33zzDdOnT+fFF1/Ezs6OypUrM2PGDBYvXlyc9YnkkdOSvGX9mjj7+HP+3Dl27NylluQiIjZkcZhITEykTp06wJVJl3Cld0dWVlbxVCZyHR4ujnwzpAVu9doBkHJ0m1qSi4jYkMVhom3btoSEhGA0Gs3bvvnmG1q3bl0shYlcT1JaBqNCw3Gt1waA1KM7GBUaTlJaho0rExEpmywOE2+99RbHjx8nICCApKQkAgIC2L17N+PGjSvO+kTyyN2SvHUTP/wbNibj4il2Rhzk8Zk7FShERGzA4kdDvby8mDFjBhcuXCAmJoYqVapoRUyxutwtyWcPa83Es3049OcBKsbuIyLKRy3JRURswOIrE8eOHaNPnz6cPXuWJk2aMHv2bPr378+pU6eKsz6RPILv9yOwibd5smVOk7lyZ/cS2MRbLclFRGzAYLJwoYhHHnmEbt268cQTT+Dg4EBWVhY//PADGzZsYPbs2cVdp8X27NlDQECArcsoE3Lmz9hyBUqTycSdd95JdHQ0MTExVKtWzWa1FLeSMN5licbbujTe1mU0Gjl48GCR/by0+MrE0aNHGT58OA4OV+6M2NvbM3ToUA4fPmzxh+3bt48+ffrQrFkzBg8eTFRU1A33f+WVVwgODrb4/FL2GAwGgoKCMJlMLF261NbliIiUSRaHCR8fH9asWZNn28aNG/H19bXoeKPRyOjRoxk+fDg7d+6kffv2NwwKa9asYdmyZZaWJ2VYzq0ONf4SEbENiydgvvbaa4wZM4aGDRtStWpVzp8/z99//83XX39t0fHbt2/H09OT3r17AzBy5Eh++OEHIiMjzetX5IiLi2PixIn069dP61jITXXq1AlPT0/WrFnD5cuXKVeunK1LEhEpUywOE+3bt2flypVs2rSJuLg4OnbsSKdOnfDy8rLo+OPHj1O7dm3za3t7e3x9ffMNE++++y5PPfUUZ8+e5fTp05aWaJZ7LQwpPunp6bYuwez+++/np59+YtmyZfTt29fW5RSLkjTeZYHG27o03tZV1ONdoEZfnp6edOjQgV69etG6dWvS0tI4c+aMRcempKTg4uKSZ5urqyupqal5tq1YsYLLly8zYMCAgpQmZVzOFa8lS5bYuBIRkbLH4isTS5Ys4d133yUlJcW8zWQyYTAYLJqE6erqSlpaWp5tqampuLu7m1/Hxsby2WefMWfOHEvLypdmA1uXrcf7VFwK61N9cXJyYuXKldjZ2eHo6Fhq25LberzLGo23dWm8b08Wh4lJkyYxbtw4evfujb29fYE/qHbt2ixcuND8Oisri6ioKGrVqmXe9scffxAbG2v+LdNoNGIymfjzzz81U1+uK2TlEX4/eomK9QI49+c2Nm/eTKv2Hc0rZQJayEpEpBhZfJsjOTn5loMEQJs2bbh48SKLFi0iPT2db7/9lpo1a+aZLxEUFMTevXvZvXs3u3fv5plnnqFXr14KEnJDOW3JM3yuPC8995cF5iChtuQiIsXP4jAxaNAgJk2aREJCwi19kIuLC1OmTGHOnDm0adOGrVu3MmnSJAACAwN1r1tuWU5b8taduwEGvp83n/CT8WpLLiJiJRavgNm5c2fOnTtnbj8OBZszYS1aAdN6StqKdTGJqdzl35z0M3/h/eSX7Jn0FN4VXG1dVpEpaeNd2mm8rUvjbV1FvQKmxXMm5s6dWyQfKFIcctqSu9VrR/qZv0g5uo1RoS10ZUJExAosvs1Ro0YNPDw8CA8PZ8WKFezYsQMXFxdq1KhRnPWJ3FTutuQBnboBYDqxi4ioBLUlFxGxAovDxMGDB+nevTuhoaEcPnyYn376iR49erB3797irE/kpnK3JV/4xgDq16/PpdPHuNsthYioBILDDti6RBGRUs3i2xwfffQRr776Kv369TNvW7BgAR9++CHz588vluJELJHTdjykX2M8XBzp168f48ePp5NDJPWadFdbchGRYlagrqE5DZVy9OnTh8jIyCIvSqQgfL3c+HpwC/PciJzVU5csCuPrwS1K1YJVIiIlkcVhwtvbm+3bt+fZtn37dnx8fIq8KJHCaNasGXXr1mXXrl2cOHHC1uWIiJR6Ft/meOGFFxg9ejRdu3alRo0aREdHs379ej7//PPirE+kwAwGAwMGDGD8+PHMnz+fV155xdYliYiUahZfmejSpQvz5s3Dx8eH+Ph46tSpQ1hYGJ07dy7O+kRuycMPPwyg+TwiIlZg0ZWJjIwMoqKi8PPzw8/vymS2RYsW4e3tXazFidyqpk2bmm91HD9+PE8PGBERKVo3vTJx/vx5HnzwQb799lvztvj4eCZPnsxDDz3EhQsXirVAkYI6FZfCmHkR9OnXH4Bff/0VuLIexei54ZyKS7nR4SIiUkA3DROffPIJbdq0YcKECeZtFStWZPXq1fj5+fHZZ58Va4EiBRWy8gjL98ew2+7KVbRffvnFvLDV8v0xhKw8YuMKRURKl5ve5ti6dSsrV67Ezi5v7rC3tyc4OJigoKBiK07kVoT0a8yZhFTCT5pwr+LD7t276T9hEX+nuKmLqIhIMbjplQmj0YibW/7P6Xt5eZmbs4iUFDldRFvcWRH7uu0B2LV+hbqIiogUk5uGibp167Jjx45839uxY4d6c0iJ5OHiyDdDWuBevwMAKUf+4JshLRQkRESKwU3DxPDhw3nttdfYtGkT2dnZAGRnZ7Np0yaCg4N57LHHir1IkYLK6SLqWLUWDhWrk372KE98sUxNv0REisFN50zcd999nD17lueffx6A8uXLk5iYiL29PWPHjqV///7FXqRIQeTuItrizor0Hf4oX3wygV1rl/F45eq61SEiUsQsWmfi0UcfpV+/fkRERJCQkEClSpVo3rw5zs7OxV2fSIHl7iI6e1hrTv1TkS8+mUD2sS1ERD1McNgBvh7cwtZlioiUGhYvp+3m5sY999xTnLWIFImru4j6+/vTtGlT9u3bx/2el9VFVESkiFm8nLbI7eLqLqIAQ4YMAeCu+HB1ERURKWIKE1ImPPLIIxgMBubNm2eeSCwiIkVDYULKBF9fXzp16sSpU6fYsmWLrcsRESlVFCakzBg8eDAAc+fOtXElIiKli8KElAmn4lLYQT0cHR2ZP38+6enpgJp/iYgUBYUJKRNCVh5h7T8pVPJrQ1xcHKtWrVLzLxGRIqIwIWVCSL/GNK/pSWatK483T/9+tnlhKzX/EhEpHIUJKRNymn+1vbcbBidXli5Zwp6jp9X8S0SkCChMSJnh4eLI1KHtcavfAVNmOsmHN6v5l4hIEVCYkDIjp/lXuSb3AZB8YA2jQsPV/EtEpJAUJqRMyN38q2279tSuUxfjmSPsjDjA4zN3KlCIiBSCwoSUCbmbf80Z3obhw4YC4H5yCxFRCQSHHbBxhSIity+FCSkTgu/3I7CJt3my5eOPP46dnR1JB9bRs2EVNf8SESkEhQkpE65u/uXj40O3bt04ezaGBytfVPMvEZFCUJiQMmvo0Cu3OmbNmmXjSkREbm8KE1JmBQUF4enpyeLFi4mLi7N1OSIity2FCSmzLqRk49OyG+np6Xmaf6lfh4hIwShMSJkVsvIIsdXbA/DdlCmYTCb16xARuQUOti5AxFZC+jXmdHwKv62qw58HD7L09w3MOuakfh0iIgWkKxNSZnm4ODJneBv8uz4EwOD/e98cJNSvQ0TEcgoTUqZ5uDjy68SXMTi5kXxkC1kpierXISJSQAoTUqYlpWXw8qK/KNeoC2RlcPnAWvXrEBEpIIUJKbNy9+to1fNhADIOriL8ZJz6dYiIFIDChJRZuft1LHrzETp27Ehy7Gl8UiLVr0NEpACsHib27dtHnz59aNasGYMHDyYqKuqafZKSknjllVdo27Yt99xzDx988AHp6enWLlVKuav7dYwcORKASlEbCGzirX4dIiIWsmqYMBqNjB49muHDh7Nz507at29PcHDwNftNmDABo9HI2rVrWbp0KQcOHGDmzJnWLFXKgKv7dfTr14+qVauyfNlSXm5XUf06REQsZNUwsX37djw9PenduzdOTk6MHDmSo0ePEhkZmWc/k8nEqFGjcHd3x8vLi169erF3715rlipl0PnkLO7s0Jfs7GwmT55s3q4VMUVEbsyqi1YdP36c2rVrm1/b29vj6+tLZGQkderUMW//4IMP8hy3ceNG/P39Lf4co9FY+GLlpkrbraePlh8ipnpH7BxmMW36dF5//XXsnNwYPieCvdGJZGdlM+lh2y1kVdrGu6TTeFuXxtu6inq8rRomUlJScHFxybPN1dWV1NTU6x4zceJE/vnnHyZOnFjc5UkZ98GDDYhJTCPB/z9c3r+aL7+dxgGvTuyNTqSZTwU+eLCBrUsUESmRrBomXF1dSUtLy7MtNTUVd3f3a/bNzMzk7bffZseOHXz//fdUrFjR4s9xdnYudK1iudIy3s7Ozsx5qg19zj7Kuv2rGf/pF9QY0YgWd1UqUStilpbxvl1ovK1L4317suqcidq1a3PixAnz66ysLKKioqhVq1ae/dLT0xk5ciR///03P//8MzVr1rRmmVKGebg48uMrD+FyV3OyLp0n5e9tWhFTROQmrBom2rRpw8WLF1m0aBHp6el8++231KxZM898CYD333+fS5cuMXv2bCpVqmTNEqWMS0rLYFRoOOVb9bnyetcirYgpInITVg0TLi4uTJkyhTlz5tCmTRu2bt3KpEmTAAgMDGTJkiUkJSXx66+/cujQIdq3b0/z5s1p3rw5Tz31lDVLlTIo94qY7Tp34W6/BhjPHGHbH1u0IqaIyA1YvQV5o0aNWLBgwTXbly9fbv768OHD1ixJBMi7IubsYa1Z6BHME088gSkijAjfRgSHHeDrwS1sXaaISImj5bRF/ufqFTEHDRpErVq1OH94Jy3d47QipojIdShMiPzP1StiOjo6mldozQoP04qYIiLXoTAhcgNdeg/AvWJVFi1axIED/zb+0qqYIiL/UpgQuYHP1x/HqUUQAO+89z7w70TN5ftjCFl5xJbliYiUCAoTIjcQ0q8x9/QaiJ2bJ2ELfmXtHzvNT3w0r+lJSD/bLa8tIlJSKEyI3ICHiyNzn+1Ew8AnwWSi95Nj8zzxocWsREQUJkRuysPFkSVfvY19+aqkHtuJ8fRhrYopIpKLwoTITSSlZfDir3/i2WEwAPGbZjPyxz1axEpE5H8UJkRuIPeqmB169qVefT+MUQfYtmm9VsUUEfkfhQmRG8i9Kuacp9oR8tGHAKT98QPhJy4SHHbgJmcQESn9FCZEbuDqVTFbdu6Bd/3mXDodyV2x282rYmrdCREpyxQmRG7g6lUxP171F7R7EgwGDi6ZSjmDUetOiEiZZ/VGXyK3s5B+jTmTkMqa8G7E7V/Nq2+MI7bxIK07ISJlmq5MiBSAh4sjs4e1psOgMRic3Jg+5Vt27NmndSdEpExTmBApIA8XR2aOvA/PewaBKZu4VV/x1aBmChIiUmYpTIgUUFJaBqNCw/Fo+SBO1epgPH2YwJFv6zFRESmzFCZECiD3uhMt7qrEkp9nY7CzZ//Cbxjw6TIFChEpkxQmRAog97oTs4e1pkfn9gwdMRpTeiqbvv+Y1xbsN++rx0VFpKxQmBApgKvXnQCwb/UwDp53kHpsBzVjdwFqUy4iZYvChEgBXL3uBMCng1rTbug4MNgx7rX/Y9veP9WmXETKFIUJkULycHFk+UdPU//+x8k0pvKfB/oTfjxWj4uKSJmhMCFSBDxcHFk9+wucvOuTHvMXiX/MU5tyESkzFCZEikBSWgbP/3KAyr1fxuDkSuK2X+j3+ld6ukNEygSFCZFCyv24aOumDZg2fQZgYseMdwj6aH6eQKEnPESkNFKYECmkqx8XHf7YIJr1HoopPYUt37zO2NlbAT3hISKll8KESCHl97jo+nnfUNW/LRkXo1jxRTBRsZf0hIeIlFoKEyKFlN/jop7uLuxasxiPO+7iwqHt+N/bn/CT8XrCQ0RKJYUJkWJS07sq69asxr5cJZIPriFh0w9MfKhJniChORQiUhooTIgUk6S0DD7ceJ6qD7+HnbM7l7b/SvtHxhKTkGp+X3MoRKQ0UJgQKQa5n/BoE9CUeQsWYefkyoX1P9DkweH8fVZzKESk9HCwdQEipdHVT3h4uDjiunQFfR4MJHbzPFoFGfH8z1Ba3FlRcyhE5LanKxMixSC/Jzx6d/8P8xevwM6lHJd2hhG77BPKO5quOVbzKETkdqMwIVIM8nvCIyktgzmRDlQbHIJ9+aqkHNrIz+88zUOfrzQvbKV5FCJyO1KYELGCvHMomhG+eweed/ljPHOEteOHcf/r04hJTNU8ClcMUr8AABZWSURBVBG5LSlMiFjB1XMomtS7i7/Ct1GjVXeykmLZ+sVYGvR6mvATF2niUyHP7ZGc2x7R8ak2/i5ERPKnMCFiBfnNoahasTyHNi6lw/C3MDg4krgllLOhr5IQfcx8XO7bHp/8fux6pxcRsSmFCREryG8OBYDBYKByi+54P/EFTnfUI/3MX2ye+BQtHhxK5JlYBk7dRkRUAk18KvDBgw3Mx2mSpoiUJAoTIjaSp9to80ZE/hlOw/7PYXBw4tjvodx9dz12LPsJU1ZmvsdpkqaIlBQKEyI2cvU8Ch+vcmz78RP+88Zs3Pw7k52cQNzv3xI7axR/LJnLEzO2cTYxLc/VityTNHW1QkRsRYtWidhI8P1+AIT0a5zn9ke5SndQpfcrpLfuT/ymH0j5Zw8pv3/Lb1tCaTK/O+Uad8PRq4Z5/1NxKXyw/BAxiWnsj04E4OvBLUhKyyA47ADB9/vh6+Vm3W9ORMoU+3feeecdWxdRlGJiYqhevbqtyygTsrKyAHBwUCa9FRVcHQls7I2zgz3w7+2LfacSaV7TkyWvBvJXuWYkVGlCeYdMLp36C2P0nySFL8PuzH5OXbzM1phs9sSk8fuh85y7ZKSJTwW+HdKC4xeSCZy8hQPRiZxPMhLY2Nv8Gf/3yz6a1KhABVetunkj+vdtXRpv68rKyuL8+fNF9vPSqn9r+/bt4+233+bEiRP4+/sTEhJCzZo18+yTnZ3N+PHjWbJkCXZ2dgwdOpRnnnnGmmWK2ER+S3DPHtaaIdkm9t9Rlxqth5B8cB2XD64h/vhBOH6Qlau+xsm7Hq51WuNSszFZ1Vpy2ZjJI9O2k5CSgZuTPW8+0CDfqxfB9/vxwfJDmIBxgf74ermZ98u9DdBVDhG5IatdmTAajQwcOJCxY8fy8ccfc/78eWbNmkX//v3z7Dd79mw2b97ML7/8Qu/evfnwww+5++678fX1tehzdGXCevSbRNFqUqMC55OMeZ76SM/MJiz8NOcuGWl45x0se38okdU6k1DRjypenlxOuEhGbBTGqAMkH1jD0TVz+WrOApLOnsQh/RJpmdlEnEtn58mE/2/vzqOiOs8Hjn+ZYUBEgRAjimIK1ogbAlIoaWMMqb9UUWlUjAset5jE+NOYxEa0VvS09ZiahpPYxGhik7Rqlp82xJ2obTUxWgU33GUdNjGswzYsM+/vD3RkiiWaYZHx+ZzDH/ed5773vc+5vDzcZW6TsxfLE86ReL6AtO8qyS2tZuxQb17bdsaqLdDHg9e2neHjo1l8c7WQ6+U1BPRy57VtZ9h5No9AHw8M1XUs+yKFYB8PYv+eQkAvdwzVdXcd4+6iI7u46p7pO8DbDYOxnt/uuNChxt1R+77TfHfEfbsX+w7wdsNYXtLxzkwcO3YMDw8Pxo4dC8C8efP4+OOPSUtLo2/fvpa4Xbt2MWfOHDw8PPDw8CAmJobPP/+c8PDwthqqEO3i5uOjjcX+PYWzOQ2XPT6YFkiXTo78bU4Y0xSczRmEd9gM6gqzqE4/iVPhZQpTz1KTc4GanAsYbvSxQ+OIztMbR48eOLp7kf7wj9jmmcfFs2XUVzii6ezOtTIj+WXV5JcZLdu+VmZk5c7zHLx4HcByw+fi/2soOG7SaTXsPpvPPy9dp6q2ocCsN5nvOmbj9BB+v/vCPdO32WRGp9Ww93xBhxp3R+37TvPdEfftXuzbbDIzb4iWltJmT3NkZGTg5+dnWdZqtfj4+JCWlmYVl56ebhXn6+tLaqp8WY+4PzX+sqsunZrW/oN7uZP85mwemzALl8il9F64hZ6z3+HB0S/TN+IZHh3+BJpOrtQV6qlOPU558k7O/30ds6c9w9dvPEfue7PJfnMCu14ZSZ+HfUlcPZvKL1Zg3LOGf7y3gs/fWknJvz6k9NtPSfvHZ6z78585tudzKlIOUnnhEKcOJTLUnIo5+zTFqaeg4BIjPEq4cCaZmvwr1FxL5fK5s/yypxFNiZ7S3HS0ZbmMexiuXr1KXUk+daXXSE/P4MS5K6RnZlFv+I56w3dkZGYxeYALuupiDIXX0FUX84y/y42YQuoNhWRk6pk8sDM6YzGGwgJ0xmImD+hMeqa+Iaa8kIysbE5eTCU9K5v68lttUwa5ojOWYCgqQFdTypRBXcjQ51BfXkR9eRFZ2bk8498JXU0phqLr6GrKmDq4a0NMRTH1FcVk6nM4fTndui07l5gAN5xqyygvLsSptoxpQ9zIyM7FVFGCqaLkjmPOXMn4Qet1xL6zcvKYPLCzXe7bvdh39rXCFp2rHJRSTV9b2AreffddMjMz+eMf/2hpmzZtGpMmTSIqKsrSNnDgQBITEy2XNY4ePcqKFSvYv3//HW0nOTmZwYMHt+zgxW3V1tYC4OTk1M4juT/czPdrCZfZe76AwN7ubJoeBMDMj0+SktdwLmLEI90orarjdE4ZnZ20VNWaMFWXU19WQA/KyMzKor6sAHNlCYMeUBQVfkdm7jXMxop22zchRFtzICnpBMOGDWuR3trsMoeLiwtGo9Gqrbq6GldXV6u2Tp06WcUZjUY6d5YbvoS4afHIHwPw+3ED6NLJkUWfp5CSZ2CItxtebs4s/eUjOGocGPXno1TVmvBw0fHlK6P430/PkpJnwL3HrcnDy9sNL8CUZ0CZ6jBXV/BjDw2mmiqu5BZirq1G1VbTw0WRW1iGMtehTPXM+ElPamvr2HIsA8z1KFM9j/f14J8Xr4Eyg1I83s8TpRSHrnwHSoEyE+zjTnJWiSUm2McNs9nMqezShjYa7h0BOJtbZhmnf4+uXLpW3rCgYEgvNwBScg2NYrpw6dqtgmiwd1cAzuUZGrW5NWl7pLsrVwpurTfoRsz5RjH9urty5XqjmJ43YvINVm1KKS7kl1v3fb3SsjywZ8OY7jbmh65nz3239/Y7et8XC2tpSW1WTPj5+fHFF19Ylk0mE3q9Hl9f3yZxmZmZ9OvXD2i4PPKfMd/H2dnZ9gGLOyb5blt9uzqzfnqIZXlZ5EA0Wo3V91XM33qyoZDorOPTuT/F080FB42DZZ0I/+4UVtRYnuwI6O3OhunDeP5vyQ1tWndCH/W3anNr9GuY0bvhD76n660/+OlOWh7qZ7IsO96I6T74VkyxkxavsFsxD96I6ZFzK6bHjbbrjdpqnbT0qL21Xs8bMd81E+N9I6awUczt2ur/Y71eN2KKGsWYnLT0bBTT+0ZMcaOYm20lzazn8wNjfuh69tx3e2/fnvpuCW12z0RYWBhFRUUkJCRQW1vL+vXr6dOnj9XNlwCRkZFs3LiRoqIi9Ho9mzdvtty0KYRo6nbv/bh5r8XXrz2Bf083y42cAb3d+Z9BXqwaN4ie7p0s8T3cO9HT3aVJWxdn6/83Ivy7E9DbnbM5ZZb+Drwy3HI5pbOTlv0vD/9BMUeXRtxTfQ/xdmPP/J92uHF31L7vNN8dcd/uxb6H3DgD11La7NFQR0dHwsLCeOutt1i7di1Go5E33ngDd3d3IiMj6dq1K/3792fIkCFkZWWxcuVKtm3bRkxMDNHR0Xe8HXk0tO3Io6Ft627y/Z9fiHXzsdP104KZOMwHdxcdgT4e5JVW07d7F34bOfC2bSt3XuCbq4UE9HZnqI8HcWMHcUpfQtp3DadSh/p4cDyzhAt5BstkVlxVi6uT9q5jpoY+zL8uX79n+g7o5UaSvoyL18o71Lg7at93mu+OuG/3Yt8Bvdz4iZem4z0aCjB48GC2b9/epH337t23BuToyJIlS1iyZElbDk0Iu3a7x059PDuzodHlktu13e4rv5dHDgSwfLHVTctHD+D3ey5a1vkhMfdS30tu3Jui0Wo61Lg7at93k++Otm/3Yt9LRv6Y0lzrpylt0WZPc7SV5OTkFrs7VTSvpqYGkHsm2orku21JvtuW5Ltt1dTUcO7cuRb7eylvDRVCCCGETaSYEEIIIYRNpJgQQgghhE2kmBBCCCGETaSYEEIIIYRNpJgQQgghhE2kmBBCCCGETaSYEEIIIYRNpJgQQgghhE2kmBBCCCGETezyDU3JycntPQQhhBDivmF37+YQQgghRNuSyxxCCCGEsIkUE0IIIYSwiRQTQgghhLCJFBNCCCGEsIkUE0IIIYSwiRQTQgghhLCJFBNCCCGEsIndFBNnzpzhV7/6FYGBgUydOhW9Xt/eQ7I7Bw4cIDIykuDgYMaPH2/5crBDhw7x1FNPERgYyAsvvEBxcXE7j9R+pKamMmTIEHJycgA5zluLXq9nxowZBAUF8dRTT3Ho0CFA8t1ajh8/ztixYwkODmbChAmcOXMGkLmkpX344YfExsZalpvLr825V3bAaDSqn/3sZ2rHjh2qpqZGrVu3Tk2ZMqW9h2VX9Hq9Cg4OVv/+97+VyWRSX375pQoNDVV5eXkqODhYHTt2TBmNRvWb3/xGLV68uL2Haxfq6urUxIkT1SOPPKKys7PlOG8lJpNJjRkzRr377rvKZDKpw4cPq8DAQFVRUSH5bgX19fUqNDRUffPNN8psNqvNmzerJ554Ql2/fl3mkhZSU1Oj4uPjlb+/v1qyZIlSSjWb35bIvV2cmTh27BgeHh6MHTsWJycn5s2bx9WrV0lLS2vvodmN/Px8Jk2aRGhoKBqNhnHjxgGQkJBASEgIYWFhODs78+qrr7J3714qKiraecQd34YNGxg2bJhlWY7z1nHq1CmMRiMvvPACGo2Gxx57jK1bt5KUlCT5bgVlZWWUlpZSX1+PUgqtVouzszP79++XuaSFxMbGcvHiRaKjoy1tzeW3JXJvF8VERkYGfn5+lmWtVouPj4/80reg0NBQlixZYlk+ffo01dXV6PV6q9w/8MADuLq6kpWV1R7DtBuXLl1i7969LFq0yNImx3nruHTpEn379mXVqlWEh4fz9NNPU1lZKfluJZ6enkRHR/Pcc88xePBgXn/9dV5//XXS09NlLmkhsbGxbNiwgYceesjS1lx+WyL3dlFMVFVV0alTJ6s2FxcXqqur22lE9i0rK4sFCxbw0ksvodFoJPctrLa2lqVLl7Jq1Sqr3Mpx3jrKyso4fPgw/v7+HDp0iGeffZb58+dTUVEh+W4FJpOJLl26sGnTJk6fPs2iRYt46aWXqKyslHy3kO7duzdpq66u/q/5be6zO2UXxYSLiwtGo9Gqrbq6GldX13Yakf06e/YskydPZtKkScyZM0dy3wreeecdQkNDrS5xgBznrUWn0+Ht7c3kyZNxcnIiMjLSMhlLvlteYmIi6enp/PznP8fJyYkZM2bg5OTE8ePHJd+tqLn5oyXmFrsoJvz8/MjMzLQsm0wm9Ho9vr6+7TcoO/T1118za9YsXnnlFRYsWACAr6+vVe6Li4spLy+nT58+7TTKji8xMZFt27YREhJCSEgIAOPGjaNbt25ynLcCX1/fJteGzWYzAwYMkHy3goKCAurq6qzadDods2fPlrmkFTU3V7fEPG4XxURYWBhFRUUkJCRQW1vL+vXr6dOnD3379m3vodmNzMxMFi5cyOrVq61u6vnFL37BiRMnOHLkCDU1NcTHxxMRESH/Tdhg3759JCcnk5SURFJSEgA7duxg5MiRcpy3gkcffRStVsuGDRswm818+eWXFBYWyrzSSsLDw0lOTmb//v2YzWa2b99OaWkpI0aMkLmkFTU3V7fIPN7Sj6S0l5SUFDV+/HgVGBiopkyZorKystp7SHZlzZo1qn///iowMNDq58SJE+rw4cNq1KhRKigoSM2dO1cVFRW193Dtys1HQ5WS47y1pKamqpiYGBUcHKzGjBmjTpw4oZSSfLeWffv2qdGjR6vg4GD1zDPPqPPnzyullMwlLeztt9+2PBqqVPP5tTX3Dkop1Xq1kBBCCCHsnV1c5hBCCCFE+5FiQgghhBA2kWJCCCGEEDaRYkIIIYQQNpFiQgghhBA2kWJCCCGEEDaRYkIIcdeqqqooLi5u72EIIe4RUkwIcR/p378/gYGBBAUFMXToUCIiInj77bcxmUwA5OXlERQU9L39xMTEcPXq1dYerhCig3Bs7wEIIdrWrl276N27N9Dw+u2XX36ZiooKli1bhre3N6dOnfrePkpLS1t7mEKIDkTOTAhxH/P39+d3v/sdW7dupbS0lJycHPr37w9AZWUlCxYsIDQ0lIiICNasWYNSil//+tfk5eUxd+5c9uzZQ1VVFcuWLSMiIoKAgADGjx/P5cuXAVi3bh1Lly5l+vTpBAUFMXXqVHJzc4GGF2fFx8cTHh5OWFgYK1asoL6+HoAjR44QFRVFSEgIM2fORK/Xf+++rF27lqeffhqTyYTZbGby5MnEx8e3UuaEEI1JMSHEfS4kJAStVsuZM2es2j/55BM0Gg1Hjhzh008/JTExkaSkJNauXYu3tzfvv/8+o0ePZtOmTRQXF7N7925OnDiBr68vGzZssPSzZ88eFi9ezLfffkvXrl1Zv349AFu2bOHgwYNs376dAwcOcPnyZTZv3kx2djYLFy4kNjaWo0eP8vjjj/Piiy9iNpub3Y+FCxdSVVXF5s2b+eijj6iqqmL+/PktnzAhRBNymUMIgbu7e5PXcHft2pXz58+zd+9ehg8fzsGDB9Fomv7/ERMTg4ODA46OjuTk5ODm5kZqaqrl87CwMIYOHQo0vLlw586dQMPbUWfOnIm3tzcAb775Jg4ODuzcuZMRI0YQHh4OwKxZs/jLX/5CSkqKpZ/bcXZ25g9/+AMvvvgiSin++te/4uTkZFtihBB3RIoJIe5zZrMZg8GAl5eXVXt0dDRlZWVs3LiR2NhYhg8fzurVq/H09LSKMxgMxMXFcfHiRfz8/HBxcbH6vHG8Vqvl5rsFi4qKrLbZq1cvAPLz8/nqq68ICQmxfFZXV0deXl6zxQTAsGHDePDBB3F2dsbf3/8usiCEsIVc5hDiPnfy5Enq6+vp16+fVXtaWhqjR49m165d7N+/n8rKSjZu3Nhk/VWrVhEYGMjRo0f55JNPGD58+B1t18vLi+vXr1uWk5KSSEhIoFu3bkyYMIGkpCTLT0JCAhEREd/b59atW9HpdJhMJj777LM7GocQwnZSTAhxH0tJSWHFihXMmDEDd3d3q8927dpFXFwcFRUVeHp6otPpLDE6nc5yWcRgMODi4oJGoyE1NZUtW7ZQV1f3vduOjIzko48+oqCgAIPBwJ/+9CfKy8sZNWoU+/bt4/Tp0yil+Oqrr4iKiqKkpKTZ/vLz84mPj2flypXExcXxxhtvcO3atR+YGSHE3ZDLHELcZ8aMGYODgwMODg54eXkxfvx4nn322SZxzz//PMuXL+fJJ5/EZDLx5JNPMmvWLACioqJ49dVXWbp0KbGxsSxfvpz33nsPb29voqKi2LJli+XJjP9m4sSJFBQUEB0dTW1tLVFRUUybNg2NRsOaNWtYsWIFOTk59OrVi3Xr1tGjRw8AgoKCeP/9960ugwDExcUxatQogoODARg5ciRxcXFWN4MKIVqHg7p5AVMIITqADz74gNDQUAICAtp7KEKIG+QyhxCiw1BKkZOTw6BBg9p7KEKIRuTMhBBCCCFsImcmhBBCCGETKSaEEEIIYRMpJoQQQghhEykmhBBCCGETKSaEEEIIYRMpJoQQQghhEykmhBBCCGETKSaEEEIIYRMpJoQQQghhk/8HED4LlfO6poUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 540x360 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "w = 5\n",
    "tend = 5\n",
    "dx = 0.1\n",
    "length = 100\n",
    "phi = 1\n",
    "dt = 0.001\n",
    "lab = Column(length, dx, tend, dt, w)\n",
    "D = 5\n",
    "lab.add_species(\n",
    "    phi,\n",
    "    'O2',\n",
    "    D,\n",
    "    0,\n",
    "    bc_top_value=1,\n",
    "    bc_top_type='dirichlet',\n",
    "    bc_bot_value=0,\n",
    "    bc_bot_type='dirichlet')\n",
    "lab.solve()\n",
    "x = np.linspace(0, lab.length, lab.length / lab.dx + 1)\n",
    "sol = 1 / 2 * (\n",
    "    special.erfc((x - lab.w * lab.tend) / 2 / np.sqrt(D * lab.tend)) +\n",
    "    np.exp(lab.w * x / D) * special.erfc(\n",
    "        (x + lab.w * lab.tend) / 2 / np.sqrt(D * lab.tend)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkkAAAGzCAYAAAA7YYPWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XlcVPX+x/HXDDAsgiCgaApS7ppL7kuXUil3xR0zl+pnlss1u92b12tWXkstb7lm2q3cSM19SU1xwxVTXBBRc8EFERVkX4ZZfn/MnRFkAIVhGfg8H48ejzjfs3w5TWfefM93Uej1ej1CCCGEECIHZWlXQAghhBCiLJKQJIQQQghhhoQkIYQQQggzJCQJIYQQQpghIUkIIYQQwgwJSUIIIYQQZkhIEkIIIYQwQ0KSEEIIIYQZEpKEEEIIIcyQkCSEEEIIYYaEJCGEEEIIMyQkCSGEEEKYISFJCCGEEMIMCUlCCCGEEGZISBJCCCGEMENCkhBCCCGEGRKShBBCCCHMsC3tCpRH586dQ6fToVKpSrsqQgghRLmhVqtRKpU0b968RK4nLUnFQKfTodfrS7saVkmn06HT6Uq7GlZH7lvhyb0rPLl3hSf3rnD0en2J3jdpSSoGKpUKnU5HkyZNSrsqVic1NRWASpUqlXJNrIvct8KTe1d4cu8KT+5d4YSHh6NUllz7jrQkCSGEEEKYISFJCCGEEMIMCUlCCCGEEGZISBJCCCGEMENCkhBCCCGEGRKShBBCCCHMkCkAhBCigtHr9RaZy804X43M9/Ps5N6BQqFAoVCUdjXyJSFJCCEqAK1WS1xcHMnJyajVaouc0/gFX5Lz1pQXcu8MVCoVLi4ueHh4YGNjU9rVyUVCkhBClHNarZZbt26RkZFh0fOW9VaAskzunYFarSYuLo7U1FR8fHzKXFAqFyFp8uTJnD59mpCQkKc+5s6dOyxYsIDQ0FASEhLw9fVl+PDhDBkypBhrKoQQJS8uLo6MjAxsbGzw8vKiUqVKFmnB0Gq1AGXui80ayL0ztKalpqYSGxtLRkYGcXFxVKtWrbSrlYPVt/N9//337Ny585mOuXPnDoGBgezcuZN27doxfPhw0tPT+eSTT5g9e7ZF63c7Po1pW8JJydAwbUs4t+PTSqzcWKbW6FBrdHmW53VuI3PHCiGsR3JyMgBeXl64urpia2uLUqmUf+SfUv3H1tYWV1dXvLy8cnxOyxKrbUnKzMzkiy++YN26dc987JdffsmDBw9YtmwZr7zyCgB//etfGTVqFMuXL6d37968+OKLRa7j7fg0ApedIDohnd0XYnmYksmBSw9Y+257vN2dirUcMJXFJGSgB/Zfum+23Ny5jfVfcvAa95IychybvXxpyDWmdG/E7N2RjPWrYyrLTq3RMWNHRJ7lQojio9frTX2QZJ0wURYZP5dqtRq9Xl+mXkVaZUvS/v376dGjB+vWrTOFnKcVExPD/v37admyZY5jHRwc+Nvf/oZery9U8DJnacg1ohPSAXiYkglAdEI6S0OuFXt59rJ9l+6z/9L9PMvNndsYwH45eSvXsdnLV5+4xatzD7L6xC0Cl53I1VI1dVM4760+bbbcSFqqhCg+2UexKZVW+cgX5Vz2z6UlRl1aklX+H7NhwwZSU1P59NNPWbp06TMdGxYWhl6vp0OHDrnKWrZsib29PSdOnLBIPad0b4Sns32ObZ7O9kzp3qjYy6f3bkLXhrnf7XZpWI3pvZsUeO4lBx+HqCePhYIDXEEhy7jPkyHqyWuChCghhBClwypD0qhRo9i3bx9vvPHGMzfLXbtm+JL28fHJVWZra0uNGjW4c+dOkYfI6vTwwbqzpgBh9DAlkw/WnSU1U1Os5VlaHebyuALI0uoKPPe9pNyjYLLf6aKGrLxC1I/Hbpv2l5YoIYQQpckq+yS1a9eu0McmJSUB4Orqara8cuXK6HQ6UlJScHd3L/R1rt68w68LFoPCkEMdVbakZxnmxfh1v4L1ixVgDHgKxf/KtYCC9SEK1i9VGEKJQgEocLK3JU2tBYWCDUdhw48KFDw+RyV7O1LVWlDAxlAFm1YoQGmLQmkDShsUNragtGFbpA3bt24BGxsUSsM2t0oOJGbqUdjYsDP2GjtDTqKws0dpZ4/C1h6FnT0KpQ37Lt1nXFAY8wNb5Buy5ge2KDBk5RWipnSrC+Tsz2VkbImaGdDUtE2t0TEuKIzgyNhcfaaEEEKIorDKkFQUaWmG1gaVSmW23Li9qC1J2vRkksN+M/1cUJ/9pALKEwsoT3iaSuXh/tPsZGOL0tae5Xb2rJjsYAhRtvYo7Cvh6OxClq0TSgdnNoY6s3ltJZQOzigdnLGpVAWbSlVQOjiz79J9xq48yZz+jYh+lJLrElqtloz0dLQ2Shbsv5UrRL1Sz52/dfElNTUVgKi4NMYEneNu4uPXfYv2XWZ6z/o5jlNrdcz+/SrvdPSmpptjoe5RWZeenvs1pXg65f3e6XQ6dDodCoUCrVZr0T4f5Wm26CVLlrBw4UIcHBw4ePAglStXLtbrPXnvpk6dypYtW/jvf/9Lx44dn/l8cXFx7Nmzh2HDhpm2jRo1ij/++INdu3ZRu3btItc5P40bN8bLy4sDBw4887E6nc40C3xqamq+fedKus9ShQtJ9vaGV0RZWVlmy43hyMlJWiNy0GrQaTWQmZqr6Kmmp7OxxaZSFX6pVIV1s92xqVQFW9dq2LpWx9bNC1u36hy8oufjbVpm9qrL/eRMMyd53BYVFZfGwGWnyNA8ftB0ru9haokyUmt1fLghggNX4jh8NY7lI1uU26AkhCgcvV7Pli1bcHR0JD09na1btzJixIjSrtZTi4uLo0ePHtSvXz9HSAoICKBNmzZ5vjkRBatwIcn4YTG+dntSUlISCoUCZ2fnIl3HxrEyLm0CQK/Hp4ojfvU8OXTlAbfjUwE9hjCsB72eWm6OdKjjzrGrD4l+ZGjpMqRlPejhOVd7WvtW4Y8b8cQkpoNej/5/ZaDHy8Wel7zdOHPrEbFJ2cp1WvRaLVUclei1Gh6lZKDXadBrtYYynRZHGz1uDjbEJaeRkan+3zEa9BpzIaUItBq0SQ/QJj3IcxeFypG1btXZ+E0t7Dx9UFWtjZ1nbWzdqqNQ2nDozzj+vvky8wNbMPaX0BwBydPZnoVvtMbZ4fFHWq3RMSkojANX4gC4m5jJyj9icryuK29kiHfhldd7p9PpTH+Z29jYFMsIN2ufEPHEiRPcvn2b999/n59++olff/2V0aNHl8i1jffO2L/Wxsbmme+nWq0mJSUlx/kABg0aZKFaPr3CfBaMa7gpFIoCJzot6ekBKlxIqlOnDgC3bt3KVabRaIiJieH5558v8oPEppIb7l3+D4AU4KSzPakumZjr5ZQOnHe2J8MzEw8z5WrgirM9WbUy8TRTrgNuOtujq5tJ1XzqVD2P7Z7O9iif6F+k1+vRa9ToszLoUNsFrTqDY5fvotdkos8y/KPLyqC+mxL/Os5sPnGZG3cfoMtMQZeRgi4jFV16Evq0BEML1FPQq9PJun+DrPs3cmy3sbPHxr0Wqhr12HK+ATtDTmLn6Y1C8fi/kbE/1HfDW6KyVXL9QQojfzzJnWyv7PwbeZk6jhvJHE5CiI0bNwLQrVs3rl69yt69ezl58iRt27Yt5ZqJ0lbhQlKbNm1QKBSEhoYyYcKEHGWnT58mMzOTVq1aFfk6z3tWwr+RLcGRsQA5Ojn7N/JifmALJq09W6zlE9ecMY0cM+rasBoLhr309Oe2s+fkAwAn7KvnfJUFcBfYaWtPcotWeLbIfR/0ej0vezuQnviQo+HX0KY+Qpf6CE1yHJrEWBwz4lA/ijH9FWSONisTbew11LHXSDm7GzC0OtnXaIB7vZboajRBVb0uwZGxjAsK4x/dGtB38REysh63NPk38jIFKCPp9C2ESElJYe/evVSrVo2GDRvSq1cv9u7dy9q1a3OFpBEjRhAZGcnevXv59ttv2b9/P4mJidSuXZthw4YxfPjwHPtnZmbyyy+/sHv3bq5du0ZGRgZubm60bt2asWPHUr9+zv6T2f3f//0fhw8fZvXq1bRp0yZHWVJSEi+//DL16tXj1VdfZdGiRYBhipsGDRowYcIEJk6cyIgRIzh58iR79uzJ0SfpwoUL/Pe//+XUqVOkpKRQs2ZN+vXrx8iRI3FwcDDtd/PmTX788UeOHz9ObGwsCoWCmjVr8vrrr/Pee+/l2Le8qnAhqXr16nTq1IkjR44QHByMv78/ABkZGXzzzTcAvPHGG0W+jlIB84a24NW5B3OEEE9ne+YNbUEle9tiLbezUWKuUVIP2NkoCzz3d8Nb8t7q00UOWQqFgqN3MgEXHH3NpCigc4OqqFMS2X/qAppHMWTF3Ub98CZZD26iSYgBMx319Op0Mm6e5e7NswAo7Cvh4P0im59vye6T7bB1edzmZvy98gpIYH7knBCi/NuxYwfp6ekMGzYMhUJBly5dqFy5Mnv27CE+Pj7XKGeNRsOIESNIS0ujW7duZGVlsWPHDmbMmIFSqTT1CdLpdIwZM4bQ0FBatmzJoEGD0Gq1hIWFsWvXLkJCQvjtt9+oUaOG2XoNHDiQw4cPs2XLllwhaefOnWRmZjJgwADq1q3LyJEjWblyJTVq1GDgwIH5toDt3buXyZMnA9ClSxdq1qxJaGgo//nPfwgLC2Px4sXY2Nhw6dIlhg8fjkajwd/fn+eee474+HiCg4NZsmQJ169fZ8GCBUW59VahXIekyMhIgoODqVmzJgMGDDBtnzZtGoGBgfz1r3+lR48eeHl5sW/fPqKionjnnXdo3LixRa4/e3ek2WHys3dHMjOgabGWg2Gm7Sftv3SfGTsiTPvmdW7AbMiqVtn+qQLcs4SsA5cN/ZTsa9THvkbOv6x0WRm0cEkl/uZlzoWdIjPmClkPb4E+58gQfWYq6VdDSb8aCnuXoKpRH6d67XGq14GHeOd4FafW6Hjr5z84eu2h6Xhzr+JAXseJiiMtLe2ZR/WWxiKtKpXKogNrjK/a+vXrBxgG9xhXdNiwYQPvvvtujv3T09OpXLky69evx9HRMAikT58+vPnmm6xevdoUkvbu3UtoaCgDBgxg1qxZOc4xefJkdu7cSXBwcJ4dxLt27Yqbmxu///4706dPNw06Ati8eTMqlYrevXvj6upKzZo1TSFp4sSJef6uKSkpTJs2DTs7O1asWEGzZs0AQ4v/+++/z4EDBwgJCaFz587MmzePlJSUXC1ZH374Ia+//jp79+4lJSWlyP13y7pyH5IWLVpE27Ztc4Sk559/nnXr1jFv3jyOHDlCZmYmvr6+zJw506Id3cb61eHApQdEJ6Tj6WzPw5RMaro5MtavTomUG8u6NqxmWrvNXLm5Y2fsiMgVsp5zdWDcq4ZXbgUFOMg/ZOUVol6p58HXAxoxddufBEfGorRz4HyGA3h1xKOHYVisTp2OOuZPMm6Hk3HzHOqYy4bO6NmoY66gjrlCQshK7Ko9z6YmXUhLGMQP7/nz7srTOQJSl4bVcr2KA3kdJyqODz74gIULF1rFkH6lUsnEiROZN29ekc919epVzp8/T/369WnYsKFpe0BAAOvWrePXX39lzJgxuToLjx492hSQwNCNw8XFhZs3b5q2NWjQgC+//NLscP727duzc+dO4uPj86ybSqWiT58+rFq1in379tGzZ08Abty4wdmzZ+nZs+czj1oLCQkhISGBUaNGmQISGDpDf/jhhzRs2NDUcjZixAi6d++eqxXLw8ODevXqcebMGRISEiQkWYPLly+b3T5gwIAc4Sg7X19fi/xPlh9vdyfWvts+z0Vgi7vcWGZsIXmyRSS/Y7MHsK4Nq1Gtsj3jXq1rtrwwIQvMh6iqzioqqZ62JaoZvDwcnTod+7g/ib0YSvqfJ9AkxuY4Juv+DR7d/5E1B39m84KXcGneHce6bQ0TbeZRD3kdJyqSRYsWWUVAAsNrrEWLFlnk+W1sRerbt2+O7S1btsTX15eoqCiOHDnCX/7ylxzlvr6+uc7l4uKSYxV7X19ffH19ycrK4uLFi0RFRXHnzh2uXLliWvqqoHs+aNAgVq1axdatW00hafPmzQB5frfl5+LFiwC0aJG760P9+vVz9JHq1KkTAAkJCVy+fJlbt25x69YtIiIiuHDhwlPVvzwoFyGpLPN2dzJ9sZr7gi3O8uxlBZWbK8sesp5sZSkooBUUsvIKUe++/Hi5mIJaoowhRqlyJKtGM9xrNEPf5f/IehBF2p8nSLtynKz71x8frNeRcf00GddPY1O5Gi4te+Lc7HX2XYJxQWE5XsdlD0iQ9+s4IcqDCRMmWE1Lko2NTa5BN4Wh0WjYtm0bAHPnzmXu3Llm91u7dm2ukJT91ZfRk61Ner2en3/+mR9//JGHDw0t105OTjRu3JhGjRrx4MGDAidGbNiwIU2aNOHIkSPExcVRpUoVtm3bhpeXlynEPIuEBMO0wy4uLgXue//+fWbPns3vv/+ORmMYoezl5UXLli3x8vLizp07ZW4x2uIgIUnk6cmQlV/5s4asvEKUu73hf7qCWqJUtuY7n3u7O/H7jPeYtLY9wZHDyIq7Q0rEflIvHECb/HiOJm3SfRIOLifxyC84N3uN3UkDGQe5Rg2C+ZFxQpQn8+bN48svv6xQfZIOHjzIw4cP8fX1zXOpq/Xr13PgwAFiY2Px8vJ6pvOvWLGCOXPm0KRJEz7//HMaNmxIzZo1USgUrFmzhpCQkKc6z8CBA5kxYwa7d++mXr16xMTE8N577xVqmhrjXGDZW7yyS0tLw8nJCb1ez7vvvktkZCTDhw+nV69e1K1b1/R6b8iQIdy5c+eZr2+NJCSJYpNfyMorRBmXHCmoJUqtyb1Ir4OtkpVvt32ipQmq+I3E7S9vknErnJQzu0i7cszU8VuvUZMc9hvJZ39n3dEu7Do+BDu3xzNK5dVfyUg6dovywsnJ6ZnDR2mEJEvZsGEDAO+//z4BAQFm97l37x6HDh1i/fr1z9x6tWXLFgDmz5+Pt7d3jrKrV68CT7fERu/evZkzZw7BwcGm+f369++fY5+nnWDR2O/q/Pnzptd3RpGRkQQEBDB06FDeeOMNIiMjadOmDdOnT8+xX1ZWFlFRUU9df2snfxqLUmMMUeYCiDFEvdnehyVvtmLWgGY5QsiMHRE5WntqVXFk56S/8EJVQydCY0uTp7OhWVyhUOJYuzlVA6ZQ8/2fcO0wFKWT2+ML6jSknN/D3WXvErd7EdqUR4bj8qm/8bXc6hO3CFx2gtvxaUW4G0KIkvLw4UMOHz6Mk5MT3bp1y3O/IUOGAIYWJe0Tg0MKYpxDKC4uLsf2kydPsn79egDTa6z8uLq64u/vz8mTJ9m1a5epv1R2tra2T3U+f39/nJ2dWb9+PZcuXTJt1+l0LFmyBAA/Pz/T68TExMQc59RqtcyaNYvExMSnrr+1k5YkUWbl1xKVvaUpr4kizbU0bZ3Qia9/v0KwiyeuHYeSEh5M5ulNpMbdM+yk15FybjepFw9Sud1A9qr7My5Ike9ElCAdu4WwJps3b0aj0dC3b98co9Se9Oqrr1K1alXu3bv3zAu3DhgwgDNnzjBmzBh69OhB5cqVuXTpEseOHaNKlSpkZmaa+ggVZODAgfz222/ExsaaHeLv7u6Ovb09Fy9eZObMmXTo0IGuXbvm2s/FxYWZM2fy0UcfMWTIEF577TWqVavG8ePHiYyMpE+fPvj7+6PX62nZsiVhYWEMGjSIDh06kJWVxeHDh4mKisLDw4O4uLinrr81k5YkYZWytzQVFGDgcUtTg+qV+W54S/wbeaGwVeHyUk883lqCR4+/Yuv2eFI3fVYGiUeCuPvDu2zd+Cvvrz6N+n9rxUnHbiGsm3GE2JOvrZ5ka2trGkW2du3aZ7rGkCFD+OKLL3juuefYvn0769evJz4+ngkTJrBr1y4cHR05fPjwU7VQdezYEQ8PDxwdHenRo0eucjs7O2bMmEHVqlVZu3YtwcHBeZ6rR48eBAUF0bFjRw4fPsyqVavIyMjgo48+Yvbs2YDh9d3ixYsZNmwYSUlJrF69mn379uHt7c2yZcv4+OOPAZ45OFojhb4ivFQsYREREeh0Opo2lVaFZ2Xsk1SUxUanbQln9YnHa/Pl1dL05BQDeq2G5LO7SDy6Bl16zgWQHZ5vRd9xn/DTxF5lsmO3Je5bRVXe751OpzNNk9KgQQOLLnBrzX2SStuz3Lu7d+/StWtX+vXrZwoy5cmzfEbDw8NRKpU0aVIyf5RKS5Iod8b61aGmm6EJPb8A82R/o780qM7tHYsY/p+tVG43CGzsTGUZN06z/p9DqdV9DHsj7pq2l4WAJIQov/R6vWl6BkssmSWejfRJEuVOQdMPGF+XZZ9i4OW6nvw0ug0qWyX/HeOHysmZXce7E//7YjKizgCg12SScHA5aZeP49nnbzz/fJ1ca8I9SUa/CSEK4+HDh7z11lskJiYSGxtL165dc8ySLUqG/PkryqX8Rs49OTLOv5GXKSCBYWTcd8Nb0qNDc6oNmYFnn7/nGAmnjrlMzM9/JfLQViatPWPqq/QkGf0mhCgsd3d3MjMzSUpKolu3buXyNZs1kJAkKpyneR1nDEpdG3lRqfErPDfmeyq9+Hi0iD4rg/jdC1gzaxLvLD2QKyjlNfpNCCGehlKpZM+ePZw9e5YFCxZQuXLl0q5ShSQhSVQ4+Y2Me5Kx35KNgzOevSbj2W8KSofHCzqmXznO2mlvMnTWWhn9JoQQ5YyEJFEh5fc6Dsz3W/Jv5EXUun8z5Mu1ONRubtquSbjH1hlv023Cl6RmaswGJOncLYQQ1kee2kI8Ia+WoO+Gt6SSvS0rJvbgjU+X4eY3EmNbk16TycGln1Cj00AZ/SaEEOWEjG4T4gnmOnZnDzoqWyVLRrRGoZzIdq86PNz+NbqMFACST28j61E0VftN4fXmvhKQhBDCisnTW4gnPEvH7j69elJ91Dzsqr1gKsu4fprYX6Zgn5VY4LXUGh3TtoTLyDchhCiDJCQJ8YSn7dhtXET3+eefp/qbX+FUv6OpTB17jaV/e4NhX22QKQKEEMJKSUgSwoyCOnbD40V07zxKR2nngGe/j3Fp3c9Urk26z5YZbzPws59ligAhhLBCEpKEKARznbtfa/IcMUc20HbYZIwdunWZqfz29UT6T/+vTBEghBBWRkKSEM+ooNFvh1fOpfO4L8HGMC5Cr05n19y/0m/q9zJFgBBCWBEZ3SbEM3qa0W+75/+D3g5OBC/4O3qNGn1WJr9/O5kXoh7i+EKrPI8VQghRdsiTWYhn9LSj33bMmcDrk79BYWcPgF6j5v6mf5N+IyzfY4UQxW/Tpk00aNCAl19+mcTEvEeiLly4kAYNGrBp06YSrF3xadCgAX5+fsV+ndDQUBo0aMBHH31U7NcqTvJ0FuIZPcvot21fvk/gJ0tQqAyhCq2GB5u/IONOJPMDW0hAEqKUPXjwgJkzZ5Z2NUrMhAkTeOutt0q7GlZDntBCFMLTjH4zcvJpiteQGSjsHADQZ2XyYMNnvPnVr3lODyCENbkdn8a0LeEFfp7L6rxg27ZtY9++faVdjRIxceJECUnPQEKSEMUkewdv+5qNqDpgmqkzty4zlW2zxzNs7uZ8v1jK6peKEEZ3HqURuOwEq0/cYlxQmNXNC9akiWFU6aeffkpCQkIp10aUNRKShCgG5kbA9enRjTVr1qJQ2gCgS0tg66z3Gblwt9kvlrL6pSJEdstCbhCdkA5AcGSs2aBUlucF+8tf/kK/fv2e+rWbsY/S+vXrc5VNmTKFBg0acOzYMdO2Ll260K9fP27fvs2kSZNo06YNL730EmPGjOHatWtotVqWLVuGv78/zZs3p0+fPmzfvj3XubVaLatWraJ///40b96cVq1aMWrUKEJCQnLt26BBA8aOHcu6devo2LEjLVq0MPUNyqtP0u7duxkxYgRt27alTZs2BAYGsnPnzlz7nThxggkTJvDyyy/z4osv0qpVK4YOHcrGjRsLvHfWSEKSEBaW3xQBgYMH0mv856bt2qQHbJ7zV9758UiOL5ay/KUiRHbTejXCv5GX6ecng5I1zAv2r3/9i6pVq7J9+/Ziee326NEjhg4dSnR0NIMHD6ZZs2YcPXqUd999l3/84x/8/PPPdOrUiX79+nHz5k0++ugjTp06ZTpeq9Uyfvx4Zs6ciVqtZvDgwfTt25dr164xZswYVqxYkeua4eHhfPHFF/j7+9OnTx9atWqVax+jr776ikmTJnH16lW6du1KQEAA0dHRTJ48mWXLlpn227BhA6NHj+bs2bN06dKF0aNH06lTJyIiIpg6dSqrV6+27I0rA2QKACEsrKApAjZ+80/8HsYTuuYbwLCEycb/fMxz7t8xZ/BLVvGlIoSRcR3D7J9ZY1CaH9iCSWvPlvl5wVxdXZkxYwbvv/8+n376Ka1atcLNzc1i54+NjaVbt27Mnz8fhcIw0ezgwYM5f/48R44cYfv27VSrVg2Apk2bMm3aNLZt20br1q0BWL16NQcOHCAgIIAvvvgCW1vDV/cHH3zAsGHDmDNnDi+//DJ16tQxXTMuLo4ZM2YwdOjQfOsWFhbGjz/+SP369fn555/x9PQEYNy4cQwYMIAFCxYwZMgQKlWqxNdff42Hhwfbtm3D3d3ddI6QkBDGjBnD5s2befPNNy1238qCsvMpFaKcKGiKAJWtkpCVX9PYf4hpW/qfJ9j+3//IZJPCKhmD0pMtSk0+/d1qPstdunShb9++xTba7f/+7/9MAQmgZcuWAAQEBJgCUvbt0dHRpm2//vorNjY2fPLJJ6aABIZw9/7776PVas1OUdC9e/cC67V161YAJk2aZApIAFWqVOFf//oXEydOJCMjA61Wy+eff86cOXNyBCSA9u3bA4YWs/JGWpKEsDDjFAFLQ64xvXcTs18IKlslp3aspkG7aG6fOwpA5J4gfPpVwaXF4wdbWf5SESI7cy1K2VnDZ3natGkcP36c7du30717d/z9/S12bl9f3xw/Ozk5AeDj45Nju729YV41tVoNQFpaGlevXsXZ2Zmff/4513nj4uIAuHjxYo7trq6uuLq6FljFHAQKAAAgAElEQVSvyMhIAF566aVcZf7+/jnugTF0RUdHc/XqVe7cucONGzc4e/YsYHgtWN5ISBKiGBinCMiPo70dZw78Rv3mbYi//ScA8Xu+w879ORx8mlnFl4oQ2alslcwPbEGTT3/PVWYN84Jlf+322WefmV53WYIxFD1JpVLle1xycjIAKSkpLFq0KM/9npwQ08HB4anqZRzR5+LiUuC+Z86c4auvviIszDAhrlKpxNfXlzZt2hAeHo5er3+qa1oTCUlClCKPKq4cP7CHxi1aoU2JB72OB1u/osaoecwP7Fbmv1SEyE6t0TFp7VmzZZPWnrWK0G987bZt2zZmzpxJ7dq1c5QbX5npdLlHpKalWX4EaqVKlQCoU6eO2dFmljp/cnIyHh4eOcrUajUKhQI7OztiYmJ4++230ev1TJkyhQ4dOuDr64uDgwOZmZmsW7fO4nUrC8r2p1WIck6t0fHVkYc551BKS+DB1llMWHVSJpsUVsPcgIPs8poeoCyaNm2aabTbwYMHc5TZ2dkBkJqamuu4mzdvWrwuzs7OeHt7c+vWLbN9fiIiIpgzZw4HDhwo1PkbNmwIwLlz53KV/fLLLzRv3pxt27axZ88e0tLSePvtt3nrrbdo2LChqbXqzz8NLeHlsSVJQpIQpSTHZJM16uP+2vuPy+5eZsN3X1jNl4qo2PIakRnxebd8pwcoq1xdXfn8c8NUHRcuXMhRZhxBdvDgwRyhYM+ePVy6dKlY6jNo0CCysrL4/PPPTX2VwNBy9emnn/LTTz+RlJRUqHMHBAQAhvmfsoewxMREVq9ejVKppGPHjqZA9PDhwxzHJyYmmjq6azSaQtWhLJPXbUKUAnNfKv0DR6Ktkcryn38CIOXMTrb8Wh/AKl5TiIpJrdExce25PEex5TU9QFn/THft2tXsxI6vvPIKNWvWJDQ0lGHDhtG6dWuuX7/OoUOHaNu2LSdPnrR4Xd555x1OnDjBrl27iIyMpFOnTiiVSoKDg4mJiaF79+707t27UOdu164dI0aMYNWqVfTt25fOnTtjZ2fH77//bhrp5+npSefOnXFzc2PdunXcu3ePBg0a8PDhQ/bv3096ejrOzs4kJyej0WhyjMCzdmX3EypEOZXfZJNLvltMq1aPO4vG/b6YnSGhef71bVy2xDjjsRAlbeZvkfkO889reoAZOyJKvK7PyvjaLTs7OztWrlxJr169uHHjBqtWrSI+Pp7vvvuObt26FUs97Ozs+OGHH/jnP/+Jk5MTGzduZMuWLXh4eDBjxgz+85//YGNjU+jzT5s2ja+++oqaNWuybds21q9fT61atVi0aBGDBw8GoFq1aqxcuZLOnTtz4cIFVq1aRVhYGH5+fmzatInu3buTlZXF0aNHLfVrlwkKfXl8iVjKIiIi0Ol0NG2a/+gmkZvxPb+xM2F5NG1LOKtP3DL9/OSXyonzl/lLx3ZoUg2jVew8fKg+6htG/qV+jhFz2cPWc672LB/Zgvo1PRHPprx/5nQ6HZcvXwYMS1IolZb721ir1XLnURpv/PcPohPS8x2Rmf3zWtPNkbXvtsfb3fyIr4rAOFy+KOGmvHiWz2h4eDhKpdK05l5xk5YkIUpYQZNN7riuxr3nh6afs+JuoT6ynLF+j2fTfbI16m5iJj8eu11Cv4EQj9WqYpgX7M32Pvm+QjO2KL3Z3qfCByRhPSQkCVHCjJNN5vWlMr13E/r06knlNv1N22JP7uD4PsPwX3Ov6zrX92BKt7ol8wsI8QTjvGAF9TFS2SqZGdBUApKwGuWnd5UQViS/ySaNf3FrNX/nl9vhqO9dBWDE6Ldp1qIlc4/F5+oD8nX/Bqhs5G8eIYSwJHmqClEGqWyVLB3Vjr4fzEFhZxh6q05LpkWXvuyNuGvaz/S6TgKSEEJYnDxZhSijVLZKgv4WwF9GfWzalnkngqQ/DAtSyrIlQghRvOTpKkQZprJV8tvCf+HUyM+0LeHwKrIe3raKtbCEEMKayRNWiDJMrdHxwbpzuL/2PjaVqhg2arN4uPMbJgadKvMzFwshhDWTkCREGZV9FJuNowvu3Sc+Lov5k80rlljFEg+idBkXZAXzi7IKUdqyfy6zf17LAglJQpRB5ob59+3ThzdHjDT9nHDkl8ezcWvly0+Yp1AoUKlUgPlFWYUobcbPpUqlKnMhSaYAEKKMyW/ZkrR+8zmwfx/R0dGg0xD327fsdf8PH2o1fDOoCeVzzmhRVC4uLsTFxREba/hMVapUySIzbxtbAMraF5s1kHtnuAepqammz6WLi0sp1yg3CUlClDEzdkTkuRaWys2N2fMWM2KwYeVudew1kk5t5YDNQGb/fpU5g18yHafW6JixI4KxfnVk8r4KzsPDg9TUVDIyMrh7927BBzwl46pWFfmLvrDk3uXk4OCAh4dHaVcjF3ndJkQZk9+yJbfj01hy1QXnZq+b9k888gvu2nje6eht2mZsjVp94haBy05wOz6tZH8JUabY2Njg4+ODh4eH6dWbJej1emT5z8KRe2egUqnw8PDAx8enTK5jZ7UtSRs3bmT16tVERUXh4OBAp06dmDx5MjVr1nyq4//44w+WLFnCuXPnUKvV+Pr6MnjwYN58802LLgApxLMyLluyNOQa03s3yTHMf2nINaIT0nF79S3Srp5El5aAXpOJ86nlPOdqWMbkydd10QnpLA25lucM36JisLGxoVq1alSrVs1iX9DlfXHg4iT3ztCKVtZb0hR6K4yyc+fO5YcffqBu3bq8+uqrxMTEsHv3bipXrsz69evx9vbO9/h9+/YxYcIE7Ozs6NmzJ5UrV+bAgQPcunWLgIAA5syZU6T6RUREoNPpaNpUvpSelTw48pc9AKVePMjD7XNNZT/99BPDR4zKsz+TzKlknnzmCk/uXeHJvSuc8PBwlEolTZo0KZHrWV1IioyMJCAggFatWrF8+XJT03FwcDDjx4+nc+fOfP/993ker9fr6dq1K7GxsaxZs4ZmzZoBkJGRQWBgIJGRkQQFBdG6detC11FCUuHJg6NgxqC09+I97q//lIwbYQB4Vq1K90/XcPh2hmlfCUgFk89c4cm9Kzy5d4VT0iHJ6p6cQUFBAEyYMCHHu3V/f3/atm3LwYMHTT3lzXn48CHR0dHUr1/fFJDA0GksIMDQGfbMmTPFVHshis64AO5rjavj/vo4FLb2ADx88IAty74y7ScBSQghisbqnp6nTp3C1tbWbEtPhw4d0Ov1nDhxIs/jK1eujI2NDffv3ycrKytHmTFcubu7W7bSQliYMSj16NAc15eHmbannN9DZvQlCUhCCGEBVvUE1Wq1REVFUb16dbMjNHx8fAC4fv16nuewt7dn0KBBPHz4kClTphAdHU1KSgqbNm1i1apV1KpVi27duhXb7yCEpahslcwPbEHl1gHYVfU1bY8P/p5vBjeVgCSEEEVkVaPbUlJS0Ov1uLq6mi03TkSVnJyc73mmT5+Ou7s733//PTt27DBtb926NfPmzcPZ2bnIddXr9TK7bSGkp6eXdhWshlqr48MNEShsbHH3H0vsmn8att+7So/3P2PHwqmobCQoFUQ+c4Un967w5N4VTkl3o7aqJ2hammGul7zm+TBuz8zMzPc8u3btYs2aNbi6ujJgwABGjBhB/fr1OXXqFP/+978LPF6I0mYMSAeuxAHg4NMUp0avmMpPrP+O8cuPyXIlQghRBFbVkmRvb+ig+mRfIiO1Wg2Ak1PeswuHh4fz97//ndq1axMUFISnpydgeJX373//mzVr1uDp6cn06dOLVFeFQiGjFopA7l3e1Bodk4LCTAEJoHN9Dz5Yv4TWrVqhyUxHl57EjhULcXL1kL5JT0k+c4Un967w5N49m5KeV8mqnpzOzs4olco8X6cZt+e3/sv69evR6/VMmjTJFJDAMNHa1KlTcXFxYePGjbJatiiT8lrX7ZtBTajr68Pnnz4O98lndrHz0HHDArga+TwLIcSzsqqQpFKp8PHx4e7du2Zbk27dugVA3bp18zxHdHQ0AHXq1Mnz/BkZGcTFxeUqF6K05bmu2//6Hv3tw8nUrVfPUKjXEb93KXsv3mPGjojSqK4QQlg1qwpJAG3atCErK4uwsLBcZcePH0ehUNCyZcs8j69atSoAN27cyFWm0WiIjo7G3t4+z87hQpSm/NZ1A8Mr6U9mPp4xPvNOBI7Rpxnrl/uPArVGx7Qt4bKumxBC5MHqQtLAgQMB+Pbbb8nIeDyzcHBwMCdPnqRLly5Ur149z+N79+4NwMKFC4mPjzdt1+v1fPvttyQkJNCtWzeLLgIphKUY13V7s72P2b5Gt+PTWHrdFccXHs8jlnl8FV7OObsfygK4QghRMJvPPvvss9KuxLOoUaMGjx49Ys+ePezevZuYmBg2bdrEwoUL8fDwYP78+aZWoNDQUDZv3kxycjIvvPACYJhLKSkpiQMHDrBp0yZu3rxJaGgo8+bNY/fu3bzwwgt8++23ODo6FrqODx48QK/X4+XlZZHfuSIxvkaVkJo3V0c7ujT0wkb5uAOj8b7NO3Cdo1fjsPN6gZSzuwA9CY8e4eHhQfv27YHc/ZqSMzRkaXV0aVgxP6/ymSs8uXeFJ/eucO7fv49CoaBatWolcj2rC0kAfn5+VKlShUuXLhESEkJcXBydO3dm7ty5ORa33bx5M4sWLcLe3h5/f/8cx9epU4c7d+5w7Ngxzp49i4ODA0OHDmXOnDlFftUmIanw5MFROMb71qXxc0TGJHMzzRZt6iPU9/4E4OTJk4wZMwYbO3uzHb+/GtQ8R+iqSOQzV3hy7wpP7l3hlHRIsroFbq2BLHBbeLLoY+Fkv2/GlqLfT10metkY9GrDpHUT/jqJ9JZvmu/4XYGnCJDPXOHJvSs8uXeFIwvcCiGKxLiuW7fWDXDtMMS0fdGixew69njxZglIQgiRvxJ5OsbExJTEZYQQ/2MMSv2Gj8GmsmFEJzoNCQeXAxKQhBDiaVjkCXn+/Hnee+89QkNDc5VpNBp69erFoEGDOHXqlCUuJ4R4CipbJYtHtqPKK6NM29KuHCMz5grzA1tIQBJCiAIU+Sm5b98+hg8fzqFDhzhz5kyu8jt37gBw4cIFRo8ezbZt24p6SSHEU1BrdExaexanRn6ovB7Pk5RwaAWT1p6VWbiFEKIARQpJ169f54MPPkCpVDJ16lRGjx6dax9fX1+OHz/OlClTAPjkk0+4efNmUS4rhChA9mH+CoUSN7+RprKMm+fYvut3Wa5ECCEKUKSQtHz5cjQaDUuWLGHEiBE4ODiY3c/e3p7Ro0czZ84cMjMzWbFiRVEuK4TIh7n13Xr37MHLf/Ez/ZwQspK9F+9JUBJCiHwUKSQdP36cNm3a0LFjx6fav1evXjRp0oSjR48W5bJCiDzktQDukjdb8dWc2Y/3i7lC+p/HCY6MlaAkhBB5KFJIunfvHo0bN36mY5o1ayaj3YQoJnkugGurpEOHDvTt29dUlhCyCr1OS3BkrCyAK4QQZhQpJFWqVAmNRvNMx6hUqiIt+SGEyFtBC+B+8cUXKBSGmbWz4m6TGnGAmm6OZhfAFUKIiq5IIem5557j8uXLz3TMxYsXZbkOIYpJQQvgvvjiiwwfPtz0c9Yf61gx6iW83Z1KuqpCCFHmFSkkdenShVOnThER8XRN9REREfzxxx+0adOmKJcVQuTD292JmQFN85wH6fPPP8fW1haApAcx7Fy/qiSrJ4QQVqNIIWnw4ME4OjoyYcIEbty4ke++169fZ8KECdjY2PDGG28U5bJCiCKwc6tOc/+Bpp9nzZpFenp6rv3UGh3TtoRzOz6tJKsnhBBlRpFCkpeXF59++ikxMTH06dOHKVOmsHPnTi5dusStW7cIDw9nx44dfPzxx/Tt25eYmBimTp1KnTrS/0GI0nA7Po3AZSe490JPbOzsAcMAjGXLluXYzzhKbvWJWwQuOyFBSQhRIdkW9QT9+vXDycmJadOmsWXLFrZu3ZprH71ej6enJ//85z/p1atXUS8phCikpSHXiE5Ix9bFA6dm3Ug+bZgBf/bs2bz77rs4OjrmmkYgOiGdpSHXmBnQtDSrLoQQJa7IIQngtddeo2PHjvz+++8cOnSIe/fukZiYSJUqVfDx8eGVV16hc+fOMqpNiFI2vXcT7iVmEhwZS+X2g0g5txu9Rm1qTXp//ESz8yxN792kFGsthBClwyIhCQzTAQwYMIABAwZY6pRCCAtT2Sr5bnjL/wUhcG7e3dSaNGv2bM44teLgtUTT/uamERBCiIpCnnxCVDDGoOTfyIvK7QehsFUBEHvvHtvWrTTtJwFJCFHRydNPiArIGJS6t2mEc/Pupu1JJzagy8qUgCSEEEhIEqLCUtkqmR/YIkdrkjb1ESnndjM/sIUEJCFEhSdPQSEqKLVGx6S1Z7F1ds/VmjR+5QlZ9FYIUeFJSBKiAnpymP+TrUnbfl3NuKAwCUpCiApNQpIQFcyTAQmge5tGvDf2XdPPSaGb2HvhjgQlIUSFJiFJiArEXEAydtL+58f/wM7ODgBt8gNSLx4iODJWgpIQosIq8ZD0xx9/sGjRIpYtW8axY8dK+vJCVGgzdkSYDUgqWyXe3t6MGDHCVJZ4YgN6nZbgyFhm7Hi6RayFEKI8KfGQFBoayuLFi7lx4wb/+Mc/WLx4cUlXQYgKa6xfHWq6GWa+NzfM/+OPP0ahUACgib9D2p8nqOnmyFg/WW9RCFHxlHhIatiwIQEBAcyaNYutW7eiUqlKugpCVFje7k6sfbc9b7b3MTsPUv369Rk8eLDpZ9vzW1gzph3e7k4lXVUhhCh1Cr1ery/tSpQ3ERER6HQ6mjaVBUGfVWpqKmBY5kY8PUvet7Nnz/LSSy+Zft69ezfdunUr8nnLKvnMFZ7cu8KTe1c44eHhKJVKmjQpmfUkpeO2ECIHD5/61Gn5sunnWbNmmd1PrdExbUs4t+PTSqpqQghRoiy2wG12Wq2WR48ekZWVRV4NVc8991xxXFoIUQS349MIXHaClEZ9IOwIAIcOHeLo0aN06tTJtF/2UXIHLj1g7bvt5ZWcEKLcsWhI0mq1zJ8/n3Xr1pGUlJTnfgqFgosXL1ry0kIIC1gaco3ohHQcajXBvlYTMu8YRrXNmjWLHTt2ALmnEYhOSGdpyDVmBsjrZSFE+WLRkLR06VKWLVuGra0t9evXp3LlypY8vRCimE3v3YR7iZkER8bi2n4w9zcYQtJvv/3GuXPnaNSkqdl5lqb3Lpn+AUIIUZIsGpK2bNmCp6cna9aswdvb25KnFkKUAJWtku+Gt2RcUBh79a2wq/YCWfevAzDzi1m49vooz3mWhBCivLHoky0mJoaePXtKQBLCihmD0muNq+PaYYhp+4YN69l1/JzpZwlIQojyzqJPt6pVq5KZmWnJUwohSoExKPXpF4CtWw3DRr2O5FNbAQlIQoiKwaJPuJ49exIcHExycrIlTyuEKAUqWyUL32hF5TYBpm0p5/egTU9mfmALCUhCiHLPon2SRo4cSWhoKIMHD2bYsGH4+vrmOaN2hw4dLHlpIYSFqTU6Jq09S6WmXUk4EoQuPQl9ViYpZ3cxaW1daUkSQpR7Fg1Jfn5+KBQK9Ho9s2fPznffyMhIS15aCGFB2Yf5K+0ccHmpF4nH1gCQdHobe9sEMC4ICUpCiHLNoiEpICDAtDimEMI6PTkPEkCvwNFsCttMRkYGutQEUiIOEGzbjXFBYRKUhBDllkVDUkGtR0KIss1cQDJ20naJGMXSpUsBSDq5GedmrxEcGStBSQhRbhXbUy0rK4srV65w7tw5bt68SVZWVnFdSghhITN2ROQ5D9Lf/vY3U0uxJv4O6df+ACA4MpYZOyJKpb5CCFGcLL52W1paGnPnzmXzZkPTvFGlSpXo3bs3f//732XVYyHKqLF+dThw6QHRCem5hvnXq1ePgIAANm/eDIDNhR1Qtx013RwZ61enNKsthBDFwqItSRkZGYwcOZJffvkFOzs7WrduTc+ePXn55ZdRKpWsXbuWUaNGSauSEGWUt7sTa99tz5vtfcy+Qvvoo49M/37v8hk6uyfK4rZCiHLLoi1JP/30ExcuXKB///5MnToVFxcXU5larWbmzJmsX7+elStX8s4771jy0kIIC/F2d8pzsdqOHTvSsWNHjh07BkDKH5vx/scbJVk9IYQoMRZtSdq5cyf16tXjiy++yBGQAFQqFZ9++il169Zl+/btlrysEKIEZW9N2rRpE9euXSvF2gghRPGxaEi6desW7du3R6k0f1obGxvat2/PzZs3LXlZIUQJ6tu3L/Xq1QNAp9OxcOHCUq6REEIUD4uGJJVKRVJSUr77JCUlYWNjY8nLCiFK0N3ETF7oPNj0848//khiYmKu/dQaHdO2hHM7Pq0kqyeEEBZj0ZDUtGlTDhw4QExMjNny6Oho9u/fz4svvmjJywohSsjt+DQCl53gYqWWqJwMr9RTUlL46aefcuxnnG9p9YlbBC47IUFJCGGVLBqSRo0aRVJSEqNHj2b37t3ExcWhVqu5ffs2mzZtYsSIEaSkpDBy5EhLXlYIUUKWhlwjOiEdpcoBh6avm7YvWLAArVYL5J6QMjohnaUh0m9JCGF9LDq67dVXX2XChAksWrSIyZMn5yrX6/WMHz+eLl26FPlaGzduZPXq1URFReHg4ECnTp2YPHkyNWvWfKrjY2NjWbx4MQcPHiQhIYEaNWrg7+/P2LFjqVy5cpHrJ0R5NL13E+4lZhIcGYtLy94kndwMeh1RUVFs27aNXn36mZ2xe3rvJqVYayGEKBybzz777DNLnrBt27b4+fmZfnZ0dOS5557j1VdfZfr06fTr16/I15g7dy5ff/01rq6u9OnTBzc3N3bt2sXWrVvp1q0brq6u+R4fFRXF0KFDCQsLo0OHDrz88svEx8eza9cujh8/TkBAALa2hc+PDx48QK/X4+XlVehzVFTGObRUKlUp18S6lNR9s1Eq6P5idSJjkolK1pP18BZZD28BcPtONCftmuc5Y3dZJZ+5wpN7V3hy7wrn/v37KBQKqlWrViLXs/iM2wDNmjWjWbNmxXFqIiMj+eGHH2jVqhXLly83fcB69uzJ+PHj+eKLL/j+++/zPceUKVOIj49n0aJFdO3aFTC0ck2dOpVNmzaxefNmAgMDi6X+Qlg7la2S74a3ZFxQGL9F9yPt0mEAjh09wvW6R7GvXhewjoAkhBD5KdLTS6fT5fr5af8prKCgIAAmTJiQI4H7+/vTtm1bDh48SGxsbF6Hc+bMGc6cOcOAAQNMAQlAoVDw3nvvERAQgJOTzB4sRH6MQamX/yuoajQwbU8+tRWQgCSEKB+K1JLUpEkTJkyYwPjx400/Pw2FQsHFixcLdc1Tp05ha2tL69atc5V16NCBkydPcuLEiTxf6x08eBCAHj165CqrXbs2c+bMKVS9hKhoVLZK5ge2YOvGvjzc/jUAqZGHcXv1LeYHdpOAJISwekUKSTVq1Mgxs3aNGjWKXKH8aLVaoqKiqFmzptn3uD4+PgBcv349z3NcunQJgDp16rBq1Sp+/fVXoqKi8PDwoHfv3owfPx5HR8fi+QWEKEfUGh2T1p7FqUEnbA7+jDb5Ieg0pJzZyaS1jaQlSQhh9YoUkvbv35/vz5aWkpKCXq/Ps2O2MbAlJyfneQ7jq7iZM2cSEhLCa6+9Rvv27Tl69Cg//PADp06dYuXKlUXuTKfX60lNTS3SOSqi9PT00q6CVSrp+6bW6vhwQwQHrsShsLHFpWVvEg4tByD5zE72nB/MWK2GbwY1QWVTtoOSfOYKT+5d4cm9Kxy9Xl+i17Po02vLli2mlpq8nD59mkWLFhXq/Glphgnp8gowxu2ZmZkFnuPIkSOsXbuWb775hn/9619s2bKFzp07c+bMmVwT4wkhHssekIx6DAg0tcDq0pNIvXiIA1fi+HBDBGpt4fsgCiFEabLo6LYpU6YwceJEGjZsmOc++/btIygoiAkTJjzz+e3t7YHHQyefpFarAfLteG1cV2748OE5+lCpVCqmTp3KgQMH2LFjB++9994z1y87hUJBpUqVinSOikzuXeEU931Ta3RMCgrLEZCMnbQnRY4yjSxNPrUV52avceBKHH/ffNkqXr3JZ67w5N4Vnty7Z6NQKEr0ekUKSVu2bOHYsWM5tu3ZsyfPBWw1Gg0hISE5+jE9C2dnZ5RKZZ6v04zb8zu/saxp06a5ynx8fKhcuTK3bt0qVP2EKO9m7IjIcx6kSZMmmUJS1sObZNw8h6NvC4IjY5mxI4KZAbn/nxNCiLKsSCGpVatWTJ8+3dSCo1AouHz5MpcvX873uMK20qhUKnx8fIiOjiYrKws7O7sc5cZwU7du3TzP8fzzz3P+/Pk8W6M0Go103BYiD2P96nDg0gOiE9JzDfNv2LAhPXr0YNeuXQDYX9oNvi2o6ebIWL86pVltIYQolCKFJG9vb3799VeSkpLQ6/WMGjWK/v37079//1z7KhQKbG1tqVGjBtWrVy/0Ndu0aUNUVBRhYWG0a9cuR9nx48dRKBS0bNkyz+Pbtm3L1q1bOXbsGH379s1RduXKFdLS0sy2MgkhwNvdibXvtmdpyDWm926S6xXaBx98YApJt88d4b0xWqYOa4+3u8w9JoSwPkXuk5S9/1H//v1NkzoWl4EDB7J+/Xq+/fZbli9fjoODAwDBwcGcPHmSrl275hvCunfvzjfffMO2bdsICAigffv2gKGz96xZswAYMmRIsdVfCGvn7e6U56uz1157jcaNG5vmQbO9tBdv975m9xVCiLJOoS/p8XQYhvAVpfPVjBkzCAoKwtfXl65duxIbG8uuXbuoUqUKa9euxdvbG4DQ0FBOnjxJo0aN8Pf3Nx1/+PBhxo8fj06nw9/fn2rVqnH48GGuX79O9+7dmT9/fpF+v4iICHQ6nbRIFYJx2gTpzPhsytJ9W7p0qfhxsLQAACAASURBVOmVurOzM9HR0WV60eiydO+sjdy7wpN7Vzjh4eEolcqnnry6qCwekh48eMDevXuJj49Hp9OZ5jTQ6/VoNBoSEhI4fPgwhw4dKvQ19Ho9QUFBrFu3jqioKNzc3GjXrh2TJk0yBSSAhQsXsmjRIvr378/s2bNznOPy5ct89913hIaGkpaWho+PD4MGDWLkyJGmEXCFJSGp8OTBUThl6b6lpqZSq1YtEhISAFiwYAETJ04s5VrlrSzdO2sj967w5N4VjlWHpCtXrvDGG2+Qmpqao7XoyX93dHTkzJkzlrpsmSMhqfDkwVE4Zem+3Y5PY/Db7xO6dSUA9evXJzIyMtcfH2qNjhk7IhjrV6dU+yyVpXtnbeTeFZ7cu8Ip6ZBk0YlLlixZQkpKCj169GDGjBl4enri5+fHZ599xhtvvIFKpcLT05PDhw9b8rJCiDLidnwagctOcLu6H/zvD6MrV66wd+/eHPupNTrGBYWx+sQtw/7xaaVRXSGEyJdFQ9Iff/xB06ZN+eabbxgyZAjt27cnLS2NwMBAPvnkE3744QcePXpEUFCQJS8rhCgjloZcIzohHTu36jjWaWPavnDhQtO/GwOScb6l6IR0loZcK/G6CiFEQSwakhISEnIMv69Xrx4RERGmn9u2bUvHjh1z/VUphCgfpvdugn8jLwBcWvYxbd+5cyfXrl3LFZDAMCHl9N4l03QuhBDPwqIhydHRMUe/g1q1apGens7du3dN2xo2bEh0dLQlLyuEKCNUtkq+G94S/0ZeOPi2wNa9FmDoizhvwUKzAckaliwRQlRMFn0y1alTh/Pnz5t+fv7559Hr9TkWvU1MTMx3AVohhHUzBqXXGlencqvepu2Ll/6XPeeiTD9LQBJClHUWfTp169aN06dPM23aNGJjY6lXrx6enp5899133L17l/Pnz7N7925q165tycsKIcoYY1DqNSAQhcqwzI8+M5XUiAOABCQhhHWw6BNq+PDhtG3blo0bN3Ls2DFsbW15++23uXDhAl27dmXo0KEkJSUxYsQIS15WCFEGqWyVLHmrE85NH0/kmhy2A71ez/zAFhKQhBBlXpGXJclOpVKxYsUK9u7dS4MGDQB4++23USqVbN26FXt7ewYOHMiAAQMseVkhRBmk1uiYtPYsLi17k3x6OwBZD2+Rces8k9ZWl5YkIUSZZ9GQFBQUxEsvvcTrr7+eY/vo0aMZPXq0JS8lhCjDso9is3OvicPzrci4cRqA5NPbCa7dnHFBYRKUhBBlmkWfTgsWLOD777+35CmFEFbG3DB/v4DHr9jTr55EkxhLcGQs44LCUGt0pVFNIYQokEVDklqtplatWpY8pRDCiuQ1D9LW2eOpW7euYYNeR3LYbwASlIQQZZpFQ1JAQADbt2/nxo0bljytEMJKzNgRYXYeJAeVLePHjzdtz4gIRpeVARiC0owdEbnOJYQQpc2ifZJq1KiBXq+nd+/e1K9fH29vbxwcHHLtp1AomDNnjiUvLYQoA8b61eHApQdEJ6TnGub/1ltvMW3aNFJTU1GnJlH70RluV+tATTdHxvrVKeWaCyFEbgq9Xq+31MkaNmz4dBdVKIiMjLTUZcuciIgIdDodTZs2Le2qWB1ZGbtwytJ9ux2fxtKQa0zv3SRXp+xx48axZMkSAJo2bUafz1fx3it18XZ3Ko2qAmXr3lkbuXeFJ/eucMLDw1EqlTRpUjJLGVm0JWnFihUo/rfytxCiYvJ2d2JmgPk/ECZMmGAKSeHh51nkkVCqAUkIIfJj0ZDUrl27AvdRq9Xcvn3bkpcVQliJxo0b07VrV/bt2wfAwoUL8fPzK+VaCSGEeRbtuN2oUSMWL16c7z6LFi1i2LBhlrysEMKKTJw40fTv/8/encdFWe0PHP/MMIysikuCu4mBRG64b9RVcruYuWSuaVpqrpmVtmjlbdF7LU0tt8pSSczUSn+uuGSpiLuEW6koIoGC7MswM8/vD2IUZ0CFYdi+79er13XOOc/znDn3mZkvZ3s2b94sfzQJIUqtIvUkXbx4kfj4eNNrRVGIiori8OHDFsvr9XoOHTpEdnZ2US4rhCjDmnfsSpWatUmKu4HBYGDZsmV89NFHZuV0eiNztkYwzt9ThuSEECWiSEHS1atXmTx5smkekkql4ueff+bnn3/O9xhFUXjyySeLclkhRBkVlZDOsK+Ponq8J8R9A8CKFSuYNWtWnpWwd++3tO/8TYLHtpdASQhhc0UKkp5++mleeuklU2/S5s2badKkCT4+PhbL29vbU6tWLYYOHVqUywohyqjlBy4RnZiBS7OnSfo9CEWfxa1bt1i/fj0jR44EzDekjE7MYPmBS/lOBhdCiOJS5Inbr7/+uunfmzdvJiAggEmTJhX1tEKIcmh2oC9/J2URcg6cfZ8i9fROABYtWswLL7xAtkGxuGP37EDbLPcVQoi7WXV12/nz5615OiFEOaPVqPlymB8Tgk6w7WYfU5B04sRx9v76G+uuOlncsVsegiuEKAlWDZJyxcfHc/36dbKzs8lvr8o2bdoUx6WFEKWcKVAC1u5+gqyoPwB4ZsIsHnnmTVM5CZCEECXNqkFSeno6s2bNYufOnRgMhgLLlucdt4UQBcsNlK4cHcHeL2YAkH7hIPqUeDSu1SVAEkKUClYNkhYvXsz//d//UaVKFVq1akXlypVlB24hhEVajZpN86ZQffV8DCk3wWgg9dQO3LoM4/PBLSRAEkKUOKsGSTt27KBu3bps2rSJypUrW/PUQohyRqc3Mv3HP3D1603ir98BkHJqO1U6DGJq8CnpSRJClDirfgPdunWL7t27S4AkhCjQ3cv8XZp1Bzt7AIzpiaRd+J2Qc7FMCDqBTm8s4ZoKISoyqwZJtWvXJjEx0ZqnFEKUM/fug2TnVIXHOvY05acc3wIggZIQosRZNUgaMGAAO3fuJCYmxpqnFUKUE/cGSJCzim3tp+/dKRNzkawbFwAJlIQQJcuqc5I6duzI7t276d+/P3369KFBgwZUqlTJYtmBAwda89JCiDJgztaIfPdB6ty5M7///jsAzn+FoK/tDeQESnO2RsiO20IIm7NqkDRw4EBUKhWKorB69WqLK9sURUGlUkmQJEQFNM7fk33nbxKdmGG2zH/y5MmmIOnasRAGPj+FQzEG6rg5Ms7fsySrLYSooKwaJE2cOFGW/Ash8lWvmhPBY9uz/MAlZgf65lm91q9fP+rUqUN0dDTZ2dk8dvsIjdr3Z5y/pzzcVghRIlRKfltii0KLiIjAaDTStKkMDzystLQ0AJydnUu4JmVLeWm3Dz/8kFmzZgFQq1YtIiMj0Wq1xXrN8tJ2JUHarvCk7QonPDwctVqNr69tnudYbJuQpKenc/LkSfbv3w9AcnJycV1KCFFOjB071hQUxcTEsGnTphKukRCiIrN6kJSYmMgbb7xB27ZtGTp0KBMmTADg+++/p2fPnpw5c8balxRClBM1a9Zk8ODBpteLFi0qwdoIISo6qwZJycnJDBkyhC1btlC3bl0aNWpkesCt0WgkMjKS0aNHc+XKFWteVghRjkyePNn078OHD3P8+PESrI0QoiKzapC0dOlSrly5wgcffMCOHTvo2fPOBnETJkzgs88+Iz09nRUrVljzskKIcsS90ePU8Wpmer148WKL5XR6I+/+FE5UQrqtqiaEqGCsGiTt3r2bLl268PzzzwOYrXTr3bs3Tz75JEePHrXmZYUQ5URUQjqDV4SS5d3dlLZu3Tri4uLylMvdlHJt6DUGrwiVQEkIUSysGiTFxsbi4+NTYBlPT0+zLzwhhABYfuAS0YkZOHl3ws65KgA6nY6VK1eayty7a3d0YgbLD1wqkfoKIco3qwZJVapU4fr16wWWuXr1qjwAVwhh0exAXwJ83FHZ2ePSopcpfenSpWRnZ+f7WJPZgbZZDiyEqFisGiS1a9eO3bt3c/78eYv5Z86cYe/evbRr186alxVClBNajZovh/kR4OOOa4teoM7Z7zY6OprgDRstBkh379othBDWZNVvlgkTJqDRaBg6dCjz58/n7NmzAOzcuZN58+YxYsQINBoN48aNs+ZlhRDlSG6g1KNNE5yadDKlvzzjQwmQhBA2ZdVvF09PT5YuXYqTkxNfffUVe/fuRVEUXn31VVatWoWDgwOLFi3Cy8vLmpcVQpQzuYFStwGjTGlZ1yPQxV4GJEASQtiGVZ/dBtC+fXv27t3L3r17CQ8PJzk5GWdnZ3x8fAgICJAt2IUQD0SrUbPuneFU++pjdDF/ApB8fAs1ek/l88EtJEASQhQ7qwdJABqNhs6dO+fZJ+n48eM4ODgUx+WEEOWQTm9kavApXP36EP9/nwGQfu5XDP96kanBp6QnSQhR7Kz+DbNt2zb8/f3ZsGGDKc1oNDJy5Ej8/f05dOiQtS8phChn7l7F5tykC2onNwAUvY7U07sIORfLhKAT6PTGEq6pEKI8s2qQdOjQIaZPn05mZiZVqlQxpev1eoYMGYJer+fll1/m2LFj1rysEKIcuXeZv0pjT7OAAab8lJP/h2I0SKAkhCh2Vg2SVqxYQeXKldm8eTP9+/c3pWu1Wt555x02btyIk5MTS5cuteZlhRDlRH77IG1e/D4aTc7sAEPyTTL+PAIggZIQolhZNUi6cOECgYGB1KtXz2J+3bp16dWrFydPnizytTZu3Ei/fv1o2bIlHTp04PXXXyc6OrpQ58rKyqJXr154e3uj1+uLXDchROHM2RphcZl/w/p1GTDgTm+S5vxO079DzsUyZ2uETesphKgYrBok6XS6+wYZWq0Wo7Fof/XNnz+ft99+G51Ox9ChQ+nQoQPbtm1jwIABREVFPfT5Pv30Uy5fvlykOgkhim6cvyd13BwB82X+kydPNpWLOX8cP5cUAOq4OTLO39P2lRVClHtWDZIaN27Mb7/9Rnq65YdNZmZm8ttvv9GoUaNCX+PcuXOsXLmSVq1asXnzZt544w0+++wzFi1axO3bt/noo48e6nxhYWGsXr260PURQlhPvWpOBI9tz/D29c1Wr3Xs2JGWLVuaXle7upfh7esTPLY99ao5lUR1hRDlnFWDpAEDBnDjxg0mT57MxYsX8+T9+eefTJkyhWvXruWZr/SwgoKCAJg0aRJardaUHhAQQNu2bdm/fz+xsbH5HZ5HamoqM2fO5PHHH8fd3b3QdRJCWE+9ak58+GxTs+X9KpWKKVOmmF5/H7SW1/zrSIAkhCg2Vg2SBg0aRI8ePTh48CB9+/bFz8+PLl264OfnxzPPPMOBAwfo1q0bw4YNK/Q1jh07hkajoXXr1mZ5HTp0QFEUQkNDH+hcH3/8MTdv3mTu3LmmSaFCiNJr8ODB1KhRA4CMjAxWrFhRwjUSQpRnVt8n6fPPP2f+/Pl06tQJJycnbt++jb29PW3atOGTTz5hyZIlqFSqQp3bYDAQGRmJh4dHnl6kXPXr1wd4oPlF+/btY+PGjUyePFkekyJEGeHg4MArr7xier148WJ0Ol0J1kgIUZ4VS/dJYGAggYGBVj9vamoqiqLk2YPpbq6urgCkpKQUeJ7bt28za9YsmjdvzpgxY6xeTyFE8YhKSCf50X+h1c5Dp9Nx48YNNmzYYNY7rdMbmbM1gnH+njIcJ4QotDI1xpQ7IdxSL9Ld6VlZWQWe5/333yc5OZnVq1djZ2dn3Ur+Q1EU0tLSiuXc5VlGRkZJV6FMqgjtFp2YwajVp7iRlEW91gFcOrQNyFnt2rdvX1MPtc5g5LUfI9h3MZ6952L59oUWphVzllSEtisu0naFJ21XOIqi2PR6Vg+S/vrrL9atW8fVq1fJzs62+IZUKhXffffdQ5+7UqVKAGRnZ1vMz+12d3LK/y/HrVu3smPHDmbOnFmkVXZCCNv6+lAUN5Jy/gBK9+4J/wRJp06d4uDBg3Tu3DlPgARwIymLrw9FMbu3DKkLIR6eVYOkw4cP8/LLL993r6TCzklycXFBrVbnO5yWm5477Hav2NhY5syZg5+fHyNHjixUHR6USqXC2dm5WK9RnknbFU55brf/9GtOfHrO40i0NRvh0KAZmVfPALBs2TL+1e1ppgadMAVIkLPX0n/6NX+gB+GW57YrbtJ2hSdt93AKGz8UllWDpC+//BKDwcDkyZPp2rUrrq6uVn1DWq2W+vXrEx0dTXZ2Nvb29nnyr127BuTs12TJwYMHSUpK4sSJE/j4+Fgs4+vrC+TsHi6EKD20GjVfDvMzPbbEtfWzpiDp559/ZviCXwiLv/OdcO9mlEII8bCsGiT98ccf9OjRg4kTJ1rztHm0adOGyMhITpw4Qbt27fLkHT58GJVKhZ+fn8VjfXx8mDRpksW87777jpSUFCZMmIBaLV+qQpRGdwdKu5XWaKrVQZ8QjaIo7Fj/DdUCxgESIAkhrMOqQZJGo6F27drWPKWZAQMGsGHDBhYsWMC3336Lg4MDACEhIYSFhdGtWzc8PDwsHuvj45NvD9LmzZtJSUlh4sSJsmeSEKXY3YHS5tZ9Sdj1JQCpZ3bj1nkY3Vt6SoAkhLAKq0YDrVq14vjx49Y8pZmWLVsybNgwgoKC6Nu3L926dSM2Npbt27dTo0YN3nrrLVPZI0eOEBYWho+PDwEBAcVaLyGE7Wg1aj4f3IJdp7uSeGA1xsxUlOxMUk7v4vN5yyVAEkJYhVW/SV577TX+/PNP/ve//5GUlGTNU+cxa9YsZs2ahVarZc2aNYSFhdG7d2+Cg4OpV6+eqVxYWBhLliwhJCSk2OoihLA9nd7I1OBTqLUOuLToZUpPOb6FyUHH0OmL9hBtIYQAUClW3HRg0qRJREZGcunSJQDc3NxMy/bzXFSlYt++fda6bKkTERGB0WikadOmJV2VMid3bylZ8fFwKlK76fRG0+RtAH3KLaKXjQGjAYAaz7xJ3/7PPfCQW0VqO2uTtis8abvCCQ8PR61WmxZZFTerDrfd22Nz+/Zta55eCFHB3RsgAfRs60tqwiDWB68DIPnoz4T4+DMh6ITMTRJCFIlVg6Tz589b83RCCGFiKUDKXcUW/sR0U5Cki7lAVvQ5QkACJSFEkcg3hxCiTJizNcJigKTVqGnVqhX+/v6mvOSjPwEQci6WOVsjbF5XIUT5UCxr3c+cOcMPP/zA2bNnycjIwM3NDS8vL/r160eLFi2K45JCiHJunL8n+87fJDoxw+I+SNOmTePAgQMAZFw8THbi3zRs+Cjj/D1LqspCiDLO6kHSihUrWLhwIUZj3tUlJ0+eZMOGDUydOpVx48ZZ+7JCiHKuXjUngse2Z/mBS8wO9DUbQuvTpw+enp5cunQJRTFS+/pegj9eSb1q+T/LUQghCmLV4bYDBw7w2WefUbNmTebOnUtISAjh4eEcPHiQTz/9lNq1a7Nw4UIOHz5szcsKISqIetWc+PDZphbnGNnZ2TF9+nTT64i9m3E0ptuyekKIcsaqQdKqVatwdXVl7dq1PPvss9StWxd7e3uqV6/Ov//9b7799ltcXFz47rvvrHlZIYQAYNSoUTzyyCMAZGRk8OWXX5ZwjYQQZZlVg6Q//viDrl27UrduXYv5devWpWvXrpw+fdqalxVCCABuZSj4BAwyvV68eDHp6ea9STq9kXd/CicqQXqahBD5s2qQlJWVReXKlQssU7lyZdMmWkIIYS1RCekMXhHK5Uc6o9HmPNPx1q1bfPvtt3nK5W4lsDb0GoNXhBKdmFECtRVClAVWDZLq1avHkSNHzCZt5zIYDISGhuZ5dIgQQljD8gOXiE7MwM7RFcem3U3pn376KXq9HjDfayk6MYOvD0WVSH2FEKWfVYOkXr16cfHiRebMmYNOp8uTl5qayuzZs/nrr7/o1atXPmcQQojCmR3oS4CPOwCV2zwLqpyvt8uXL7Nx48Z8N6Oc2aNxidRXCFH6WXULgJdffpmQkBCCg4PZvn07zZo1w9XVldjYWC5evEhKSgre3t6MGTPGmpcVQgi0GjVfDvP7JxACZx9/0s7uB2DuvHnsymjEnvNxpvK5ey1lZ8lwmxDCMqv2JFWqVIm1a9fy/PPPo9Pp+O2339i2bRvHjx9Hr9czaNAggoKCcHR0tOZlhRACuBMoBfi4U7ldf1P6qZMn2bpjl+m1pc0ohRDiXlbfTNLFxYUPPviAd999l8jISFJSUnBxcaFRo0ZoNMWywbcQQpiYepSAoP1+ZF45AUDykY04NmwhAZIQ4oEV27eEvb09jz32GH5+fjRo0EACJCGEzWg1aj4f3ILK7QaY0jIjT6KLvczng1tIgCSEeCBW+aa4ceMG//nPfwgPDzfLUxSFrl27MnXqVCIjI61xOSGEKJBOb2Rq8Ckc6jdD63FnYnbSkR+ZGnwKnd7yClwhhLhbkYOk06dP07dvX77//ntCQ0PN8q9cuUJ8fDw7d+6kf//+8kgSIUSxunsVm0qlonK7gaa89PO/s/3gCSYEnZBASQhxX0UKkmJjYxkzZgxpaWmMHDmSZ555xqxMo0aN2LdvHyNGjCA9PZ2pU6dy8+bNolxWCCEssrTMv0/fZ3nMyzvnhWIkKXQDIediJVASQtxXkYKkb775htTUVD7++GNmzpyJu7u7xXK1atXinXfeYebMmSQnJ8uz24QQVpffPkjLRrThnbffMqWlRexDnxR3J1AySKAkhLCsSEHSb7/9RtOmTXn22WcfqPzIkSNp1KgR+/fvL8plhRDCzJytEWYBUu4qtqFDh9KwYcOcDKOBpCMbAQg5F8vcnX+VQG2FEGVBkYKk6Ohomjdv/sDlVSoVrVq14vr160W5rBBCmBnn70kdt5w92O5d5m9vb8+MGTNMZTPCd6NPiaeOmyNjOspjkoQQlhUpSNJqtdjZ2T3UMa6urrIdgBDC6upVcyJ4bHuGt69vcR+kUaNGUbt2bQAMeh11r+8heGx7U2AlhBD3KlKQVKtWLS5fvvxQx/z111/5zl0SQoiiqFfNiQ+fbWpxHyQHBwfeeOMN0+vwkB9xNKbbsnpCiDKmSEFS586dOXz4MFFRD/YU7aioKA4ePEizZs2KclkhhCiUXgOG4VS5KgDp6eksXLjQYjmd3si7P4UTlSBBlBAVWZGCpEGDBgEwadIkEhMTCyx7+/ZtJk6ciNFoZMiQIUW5rBBCPLSohHRGB4WjbdHHlLZ48WKz767cVXJrQ68xeEWoBEpCVGBFCpIaNmzItGnTuHDhAj169GDx4sWcPn2alJQUDAYDCQkJnDx5kkWLFtGzZ08uXrzImDFjpCdJCGFzyw9cIjoxA1e/QNSVnAFITk5m+fLlpjL3biMQnZjB8gOXSqS+QoiSp1IURSnqSb799ls+/fRT9Hq9xXxFUXB0dOSVV15h7NixRb1cqRcREYHRaKRp06YlXZUyJy0tDQBnZ+cSrknZIu12f3cHQIm/rSXpUDAAVatW5ezZs1SrUdPiPkvyMNz8yX1XeNJ2hRMeHo5arcbX19cm17N7//333y/qSVq0aMEzzzyDq6sr6enpGI1GsrOzqV69Oj4+PgwZMoQ5c+bw5JNPWqHKpd/NmzdRFEUmqBdCdnY2kLNyUjw4abf7s1Or6PmEB+diUohSPULKif8Do57MzEwqOTnzw/XKEiA9JLnvCk/arnDi4uJQqVTUrFnTJtezSk+SyEt6kgpP/roqHGm3B5fbo7Rh2TySQ38EQO3gQp3xX5uG4SRAejBy3xWetF3h2LonSb4BhBAVilaj5sthfvQZPg6VNmePJGNmKsnHfgEkQBJC3CHfAkKICkerUbP8paeo3LqvKS356E8YMlP5fHALCZCEEIAESUKICkinNzI1+BSubZ5F9c8Qm5KVRkrYZqYGn0Knl4feCiEkSBJCVDB3r3Kzc3Chcps7D+hOPv4LO0/8yYSgExIoCSEkSBJCVBz37oME8O/BL+Lm5gaAossg+cgmQs7FSqAkhJAgSQhRMVgKkAJ83Fkyoi2vvvqqKS3lxBYMaYkSKAkhJEgSQlQMc7ZGWN4HyU7N+PHjqVGjBgBKdhZJR3K2Bgg5F8ucrRElUl8hRMmTIEkIUSGM8/ekjlvOkv97l/m7uLgwY8YMU9n0k9vQJ9+kjpsj4/w9S6S+QoiSJ0GSEKJCqFfNieCx7Rnevr7FfZAmTJhArVq1ADDoddS4+AvBY9tTr5pTSVRXCFEKSJAkhKgw6lVz4sNnm1rcB8nJyYkpb7xteh2+72eSY65YPI9Ob+Tdn8KJSkgvtroKIUqeBElCCAFEJaTzS7o3mmp1ADAajbz99ttm5XIngK8NvcbgFaESKAlRjkmQJIQQwPIDl7iRoqOq/0hT2i+//MLvv/9uen3vCrnoxAyWH7hk87oKIWxDgiQhhABmB/oS4OOOo1cHtLW8TekzZsxEUZR8txCYHWibB20KIWxPgiQhhODOg2+fftyDqk+NMqUfOnSQDZt+shggyYNwhSjf5NMthBD/yA2UAnsE4OjZxpQ+fNyr7I64YXotAZIQFYN8woUQ4i65gVLPUa8BKgCy46+RGr4HkABJiIpEPuVCCHEPrUbNmjcG4vxEV1Na0m9rMGal8/ngFhIgCVFByCddCCHuodMbmRp8Crcuw1FpKgFgSLtNUugGpgafkue5CVFBSJAkhBB3uXsVm6byI1Ru29+Ul3z0J7YfPi0PvhWigiizQdLGjRvp168fLVu2pEOHDrz++utER0c/8PE7d+5kxIgR+Pn58cQTT9CjRw8WLFhAerpsDCdERWVpmX/fka9Qq3btnBeGbBL3f0vIuVgJlISoAMpkkDR//nzefvttdDodQ4cOpUOHDmzbto0BAwYQFRV13+OXLFnClClTuHjxIj179mTo0KE4ODiwbNkyhg8fLoGSEBVQfvsgrRzdiXlz55rS0i/8TmbUHxIoCVEB2L3//vvvl3QlHsa5c+d48803adWqFcHBwfj7+9OjRw+8vb3ZuHEjUVFRBAYG5nv8lStXmDZt+uWixwAAIABJREFUGh4eHvz888/06dOHLl268PzzzxMTE8OBAwfQaDS0a9eu0HW8efMmiqLg7u5e6HNUVNnZ2QBotdoSrknZIu1WeLlt99GOi/x0yvIy/6ZNm7J9+3Zu3MjJz467gkuz7lyJT+d2uo6uTSrmZ13uu8KTtiucuLg4VCoVNWvWtMn1ylxPUlBQEACTJk3Kc3MFBATQtm1b9u/fT2xsbH6Hs2vXLgwGA6NHj6ZGjRqmdLVazZQpUwDYt29fMdVeCFFajfP3pI6bI2C+zF+tVrNw4UJTWV3sJVLD91DHzZFx/p4lUl8hRPErc0HSsWPH0Gg0tG7d2iyvQ4cOKIpCaGhovsf7+fkxZcoUOnXqZJaXG3SlpaVZr8JCiDKhXjUngse2Z3j7+hb3QerYsSPP9B9oep15aA3LnvOiXjUns3Pp9Ebe/SlcHn4rRBlXpoIkg8FAZGQkHh4eFrso69evD8Dly5fzPUebNm2YOHEinp7mf/3t2LEDAG9vb7M8IUT5V6+aEx8+29TiPkhRCelEPzYAlb0DAOnJt/nifx+alcud27Q29BqDV4RKoCREGVamgqTU1FQURaFKlSoW811dXQFISUl56HPfuHGDxYsXAzB06NDCV1IIUS4tP3CJW7hSpdNgU9qKFSs4evSo6fW9k7+jEzNYfuCSzesqhLAOTUlX4GHkrjrLb6JbbnpWVtZDnTcuLo7Ro0dz+/Ztnn/+eTp06FC0igKKosiwXSFkZGSUdBXKJGm3wnvQtpvetSHRCWnsNfQlLXwP2fFRKIrC2LHj+PXX/RhQ8dqPEey7GG865l9e1ZnetWG5/S6Q+67wpO0KR1EUm16vTPUkVaqUs/Nt7qqAe+l0OgCcnMznCOTn8uXLDBkyhCtXrtCtWzdmzZpV9IoKIcodrZ2azwb60tXHg2rdXzGlnzp1kmUrv7IYIH020BetXZn6mhVC3KVM9SS5uLigVqvzHU7LTc8ddruf0NBQpkyZQlJSEn369OGTTz7B3t7eKnVVqVQ4Oztb5VwVkbRd4Ui7Fd6DtJ0zsPyFtkyw0xB8ehdpZ/cD8NY7s6n98nLsnN2AivcQXLnvCk/a7uGoVCqbXq9MfYK1Wi3169fnxo0bFnuTrl27BkDjxo3ve64tW7bw0ksvkZSUxEsvvcT//vc/qwVIQojyS6tR8+UwP54ZOwOVNqfX2piVxu29XwEVL0ASojwrc5/iNm3akJ2dzYkTJ8zyDh8+jEqlws/Pr8BzbNu2jTfffBODwcD777/PG2+8YfPoVAhRdmk1alaOD8DNf4QpLe3sftIvHeXzwS0kQBKinChzn+QBAwYAsGDBAjIzM03pISEhhIWF0bVrVzw8PPI9/sqVK7z99tsYjUY++eQThgwZUux1FkKULzq9kanBp3Bt2RttrTtbhiTs/JJXVh2UR5UIUU6UqTlJAC1btmTYsGEEBQXRt29funXrRmxsLNu3b6dGjRq89dZbprJHjhwhLCwMHx8fAgICgJzntmVkZFCnTh2ioqJMy/7v5uLiwosvvmiz9ySEKDvuXuavUttRvdcUYr6dCkY9hpSb/LTyf2gdnWXITYhyoMwFSQCzZs2iUaNGrF+/njVr1uDm5kbv3r2ZOnUq9erVM5ULCwtjyZIl9OvXzxQkHTlyBIDo6GiWLFli8fzu7u4SJAkhzFh6CG5v/7a4136bjz+cA0DKiW1s3f0kgARKQpRxKsXWmw5UABERERiNRpo2bVrSVSlzcveTkRUfD0farfAetO0sBUi5k7Qx6mnp58fZiAgANNXqUvvFRTzdtF65DpTkvis8abvCCQ8PR61W4+vra5Prlc9PrhBCWNmcrREWAyStRo1Wq+Xrr74yLQDRJ1wn8fcgQs7FMmdrhNm55NluQpQNEiQJIcQDGOfvSR03R8DyMv86Xs3w6Njf9Dr5yCZcb//JOP+8z4mUZ7sJUXZIkCSEEA+gXjUngse2Z3j7+haH0JYfuISm3RA01er+k6KQtONz3OwNpjLybDchyhYJkoQQ4gHVq+bEh882tTjHaHagL92bNaBG4GugysmPuhbJ1KmvAvnPaZodaJu5FUKIh1cmV7cJIURpk7sT9wTgx0vPk3RwHQCrVn1Dz38Hsiu1Xr5zmoQQpZN8OoUQwkpyA6X+L05B6/GYKX3IC6PZeeyC6bUESEKUDfIJFUIIK9Jq1Cwb2ZY+Uz5GpdECYExP5NbWT1EUowRIQpQh8ikVQggr02rUfPfas7g9dWdT2szIkyQf2SjPdhOiDJFPqhBCWJnp2W5+gTg+1t6UnnhgDc/P+Vae7SZEGSFBkhBCWFGeZ7upVFTv/Sp2lWvmZCpGdix+i9HL90qgJEQZIEGSEEJYiaVl/j1aerJ7y0ZUdnYAGFJusWnB27yy5li+gZLsyC1E6SBBkhBCWEFBz3b7l39nZrz7gSk949JRNq1awoSgE2aBkuzILUTpIUGSEEJYQUHPdgNQNwvE0bONKT/xt7X8smVLnme7yY7cQpQuEiQJIYQV3O/Zbu8905QBr81F41brnxSF+K3zeco9G5AduYUojSRIEkIIK7jfs920GjVfj32KZ95YgMreAQBjVjoD+/fnRlx8vkN1sl2AECVHPn1CCGElBT3bDXICpXVvPkfXcXfmJyX9fZXGnXqzO+KGKU0CJCFKB/kECiGEDWk1arYteJ0Wz4wxpWX8Fcbtfd8AEiAJUZrIp1AIIWxMq1FzIPhLHBu3M6WlHPuZ5KM/y47cQpQi8kkUQggb0+mNTPvhDDX6vI7Wo7Ep/fber2g95iNSM/X3PV72URKi+EmQJIQQNnT3Kja11pGaA967syM3ChfWf0zrKV/mGyjJPkpC2I4ESUIIYSMWd+Ru04QjB0Kwc3ABQNHr+HPtLNq9scosUJJ9lISwLQmShBDCBgrakbtV86Zs3fILKjt7AIwZyZxf9Sbt31prCpRkHyUhbE+CJCGEsIH77cjdM+BffPfdalDlvDamJXL+6zd56v0NpGXpZR8lIUqApqQrIIQQFcE4f0/2nb9JdGJGvgHOiGGDSc9IZ/zLOdsDGFJucnr5a3inZKNxrW4qJwGSELYhnzAhhLCB++3InWvcS6P5bOEi02t94t/EBr+DPuUWIAGSELYknzIhhLCR++3InWva1Ml89Mk802t9wnVig2aQnfg3/xvYrMDjZXsAIaxHgiQhhChldHojl2t1xa3LCFOaPimW2KAZdHl7rWwPIISNSJAkhBClyN2r2Kp0fJ6q/7rz+BJDajznvnqNttO/ku0BhLABCZKEEKKUsLTMv//I8Xy26AtQqQAwpidx/qvp+I2fL9sDCFHMZHWbEEKUAgXto6TVtKaKqytjxowCoxFFl86fq9/BNzGOU0FzeePHM7I9gBDFQD5BQghRCtxvH6XRo0bww4ZNqO0dcgooRq79vICG3Yax+2xMvscJIQpPPkVCCFEKjPP3pI6bI5B/oPNc/77s378fe9dqprTkIz9y86dPMGalU9fNkYXPt5DVb0JYiQRJQghRCjzoPkpdOrbj+NEw7B9paErLuHiYmO9e5fKf53h1/Sl0eqPFY2X1mxAPR4IkIYQoJR5kHyWd3sjnobfxGPZfHBq1MqXrb9/g7zXT+XnjeiYEnTALlGT1mxAPT4IkIYQoI+4OdNSVnKg5YDZVOg4x5SvZWdzaMp/1C2bx8jcHTYGSrH4TonAkSBJCiDLAUqDztG9toveupvu0hagdXEzpqad38P3M5xn4wXfycFwhikA+IUIIUcoVtD2AcyUNW/47mYH/WYvWo7EpX387hi0fvUSdriPY/Ue0Kb1rk5oFBkgysVuIOyRIEkKIUu5+2wNoNWrWvNqHoR+upnKHQaD656tdMZJ0KJiY714l8/o5AFQFXEcmdguRl93777//fklXory5efMmiqLg7u5e0lUpc7KzswHQarUlXJOyRdqt8MpC2z1W05VdEbGkZOrzHSqzU6vo3awOkdpG3HBuTFZUOMbMVACM6Ymkhe/GkJpAdKX6XLiVRc8nPLBT3wmZ7u2tSsnUk20w0rVJ/t9jZaHtSitpu8KJi4tDpVJRs2ZNm1xPgqRiIEFS4ckXR+FIuxVeWWi7Ko72dH/cnWyDkf8ObJ7vUJmdWkXXJjXZdiUbY+OnMGamoPv7L1O+LvYvUsNDiExW+EtfnV5Na2OnVqHTG3lx1VF+vXjTVDbAx53/DmyeJ5CCnGDqvV/+4LGarjj+88yG0tx2pVVZuO9KI1sHSTLcJoQQZcCDbg/w6vpTXL+dgbqSE9V7TMJ92H+xr9HAVMaYnkTCziWsnfE8gW8sIjUzmxdXHeXgpVumMvnNW7p3OC46McP6b1SIUkSCJCGEKAfym9x96etpDJ8bjNtTo1BpKpnysm9dY/fCadTwbs2eXw/kOZeleUuW9ln6+lBUsbwXIUoLCZKEEKKMu9/qt2Uj2zJg1ERqv7wM5ye6cXcYlHXtDLFBb/L3urfIiDyFoijsOR+XZ0PK3OG4e88/s8ed1XR310VWx4nyQuYkFQOZk1R4Mk5fONJuhVce2u69X/7gp1M3TK/vndxtp1bR8wkPLiQY+LtqU5wea48+8W/0iX+bjjEkxZEWsZeMKydQa52IMlThfGwa3Xxq8tJ3x8yG45YOb4XKaADutF1usPbTqRvsioil++PuVHG0t0UTlDnl4b4rCTJxuxyQIKnw5IujcKTdCq88tN2Drn7r+YQHETeSuZapxeWJrlSq+zj6pFgMyXcmbBtS4km/cJDUM7v460Y8K89kcj1VyXOuRjWcCWxWG6NBD+S03eWbqTz7xUFOXU8E8l8dd/fE74ocQJWH+64kSJBUDkiQVHjyxVE40m6FVx7a7kFXvxmMCv93JoYrt9IAsHfzwKXZ0zg0aIYhJT5Pz5KiyyDrWjjJx7agi7uESqNF4+aBSq3myq00zsWk0M27OnZqFdeTdPRe9BsJ6dmm4y2tjpOepjvKw31XEiRIKgckSCo8+eIoHGm3wisvbVfF0Z6uTdzNluznyg1Q9pyPM6UF+LizZXJnLmc4EefeFkfP1ih6HYaE6yjKPw/IVYzo46+Tfu4AWeE70SXdRKWpxNUsBy7GpfOkV3UGrzyWJ0Cq4VKJ4LHtcdTamdKkpymv8nLf2ZqtgySVoijK/YuJhxEREYHRaKRp06YlXZUyJy0t5y9cZ2fnEq5J2SLtVngVoe0Kmtit1ajN8g1pt0k5tYPUU9sxpCZYPKfaqQpOjdvh5NUBhwbNUWny/tjfff7LN1Ppveg3MrONFvMBohLSWbr/En8nZ7L3fBx13BwJHtueetWczN7LnK0RjPP3NMsrSyrCfVccwsPDUavV+Pra5uHM0pNUDKQnqfDkr6vCkXYrvIrQdg86sTviRjJXbqWh1jriUL8prq2foVLtJihGA9m3Y0C5E+Qo2VnoYi+RdvZXko/+RGbUHzjoU9Gr7FE7V+HKrfScITmfmgz48lCBPU1RCekMXhHKoUvxpqHAu3uaohLS+d/O87RtWJ3J606aDdfdnT9na4RZL9T98nPZsherItx3xUE2k3xAGzdupF+/frRs2ZIOHTrw+uuvEx0dff8D/3HhwgUmTpxIp06daNmyJc8//zy7du0qxhoLIUTJGOfvSR03R8A8QLrbvQN1Xbzcubz2HYbM/Ix6k9ZQrcckHB71Q22nyVNO0WeRGXmS6J0riPl2CteXjCBu04ds/PYLPEd/yrWbiXnK30rN4tX1p0xbDCzdf8lsY8quTWoyO9DXFECtDb3GU/P359mnafmBS2b59z537n75uWXe3hTO+LXH881/96dwUjP1Frc3KM783Dyd3mhxe4WSrFtJvHe90baDX2VyuG3+/PmsXLmSxo0b89RTTxETE8OOHTuoXLkyGzZsoF69egUeHx4ezgsvvABAnz59cHBwYPv27cTFxfHuu+8yYsSIItVPhtsKT7qgC0farfAqSttFJaSz/MAlZgf65ruT9t3DcZ0b1+CbUW0sDscZM1NJv3SUjIuHybhyHCU7q+CLq+3Q1nwU1zqNMbrVR1vzUexrPkqPlp58PrgFk9edZO9dc6UAuv2zzcCcrRGsDb1mdsrcYC+//OHt6/Phs01596fwAvNzg6h7gzRL+TVcKnErNSvPUGBh81eNaEYdN0cSslT5Hg+Y8ro1qYkCeYYi7863Zt1KQ35+731Jr0eoW9WRFs1s8/ta5obbzp07x5tvvkmrVq0IDg7G39+fHj164O3tzcaNG4mKiiIwMLDAc4wbN45bt27xww8/8Nxzz9GlSxeeffZZdu/eza5du+jfvz8uLi6FrqMMtxWedEEXjrRb4VWUtitoYrel4bgVL7TOdzhOpdGirfkozj5dqNymPw4Nm2PnWgPFoMeYdhu4529vRcGQmkBGzF9kXjlO2h97SD7yI2f2/MSn36zn/Jnj6BOiMWQkgVFBrXXkSkIm52JSeK+PLxtPRJOuM5hOd/dwXduG1fnh2HWz/K9eyAnw7pc/b/t5Dl2Kz1Pdrk1q8r9/Vub9b+d5Dv6Vk597jruHAgubrzMoPPlYdRbuu5zv8cev3jblXbmVZjYUeXe+NetWGvLze+/dGzmioOLRurWwhTI33BYUFATApEmT8nypBQQE0LZtW/bv309sbGx+h3P8+HHOnTtHz549adKkiSm9WrVqTJgwgaysLDZv3lx8b0AIIUqZwg7HdXjUjaPv/IvAHk9T1X8EtV74lDpTvsdz2Bwqt3+OSvWborKvZHaeXIaUm2RGniT15DZu7/2Kmxv/w42vX+HapwO4/sVI1rw1jDp+T3Fx8+ckhW4g9Y+9ZFw9TUzkX4xbuZektExeXX+KW6l5e7Jyh/PSsvT3zf87ObPA9zmzpw81XPK+hxoulZjZ06dI+dMDGt33+NmBvnRrYj73JncosrjqVhry83vvTlo7HnHJ/56yNs39i5Qux44dQ6PR0Lp1a7O8Dh06EBYWRmhoKH379rV4/PHjx01lLR0PEBoayvjx461YayGEKL3qVXMieGz7+w7H3b19QOfGNVg86HG0dmq+HOZnGo6zc3BBX9ePqnX9AHImfd+MJOvvS2TfvIJzajQZsZdJSsw7TykvBUNqPIbUeIixXOJ74PspoKrkjJ2jK2pHVxxc3DDYO6N2cGbDr45s/NoRldYBtdYRlb0jblUqk2KwQ611YntCNDuO/IHKzh6Vxh6VnT3Y2aNSqUyPZfl8cIsCg6yi5M/YfI55/Xx4+5dzBR5vaT6MCsg2GIutbqUl3/JcIBVKPjnFoUzNSTIYDPj6+lKnTh327Nljlr9161amT5/O+PHjmTZtmsVzzJgxg59++om1a9fSpk2bPHmKotC8eXPc3Nw4cOCAxeMfhMxJKryKMj/E2qTdCk/a7v7undeT29uUnZUzj8fZ2ZnUTD1Pzd+f50evblVHdr7qz9TgU3nmO1V31hIbE032zUiyb8egv32D7NsxaNPjSLt1A4PhztCYzdlpcgKn3P802n8CKA32Gjv0ihqVWg0qNajtUKnUoFajUtuhtdegM4JKZZeT9k8ZVDnHOGg1OdsgqHL7qlT//FMFKhWO9nZk5MkHVCpU/+Tnvs4pbyoAKhVOWg0ZOoPlfMBJqyFdp8/3befJV5kPyTpp7fIMWT5svrNWQ5ql6/9zrXzz7/H1jBHYOVbm2a7t71vWGspUT1JqaiqKolClShWL+a6urgCkpKTke47k5GQAi+dQqVS4uLgUePyD0Ol0GI1GwsPDi3Seiig3ZldZ+JCK/Em7FZ603f0NeFShfdVH0BuMOFXS4FG5EhfORZjaDpWKv5OzeLejC+DyT5KKelUduHzxHFP8KjGosTvpWXf/CDb557+8nLR2GA160jN1YDSgGPX//K8BtWJEpRjI1uuh7Px9L6yodjVXDCq7+xe0kjIVJKWn5ywNzG+CZW56Vlb+Ky0e5ByJBXYD359arc7zv0IIUZZp1Tm9QrfTs3nEpZJZR0NcSlaeAEhjp6ZOVUe0dne+A2tVcSDyVjoGozFPuQbVnfg7KZO0f47P6Y1Qoc5nLpOdWo3qrnPkcq6kwaOKQ55zPWx+rSqOKChciks1y/es6YJapcKolFw+UGrrZqv3blDZkWmwXYBcpoKkSpVyPjS5q1HupdPpAHByyn8X1gc5R0HHP4jmzZsX6XghhChLHnTv4/wmIFhzYsL9zvUg1yrqOYozvzTXrbjzS2ICS5nq6nBxcUGtVuc7HJabnjvsZknuMFvusNvdFEUhNTW1wOOFEEIIUTGUqSBJq9VSv359bty4YbEn6Nq1nImFjRs3zvccnp6eAERFRZnlxcbGkpWVZSojhBBCiIqrTAVJAG3atCE7O5sTJ06Y5R0+fBiVSoWfn1+Bx0POMv97HTp0CIBWrVpZqbZCCCGEKKvKXJA0YMAAABYsWEBm5p1NwEJCQggLC6Nr1654eHjke7yfnx+NGjVi69atnDlzxpSekJDA0qVLqVSpEgMHDiy+NyCEEEKIMqFM7ZOUa86cOQQFBdGwYUO6detGbGws27dvp2rVqgQHB5ue3XbkyBHCwsLw8fEhICDAdPyxY8cYPXo0KpWKwMBAXFxc2LZtG3FxccyePZthw4aV1FsTQgghRClRJoMkRVEICgpi/fr1REZG4ubmRrt27Zg6dWqeh9suXryYJUuW0K9fP+bOnZvnHOHh4SxatMg0bPfYY48xZswYnn76aZu+FyGEEEKUTmUySBJCCCGEKG5lbk6SEEIIIYQtSJAkhBBCCGGBBElCCCGEEBZIkCSEEEIIYYEESUIIIYQQFkiQJIQQQghhgaakK1DebNy4kbVr1xIZGYmDgwOdOnVi2rRp1KlTp6SrVupMmzaN48ePc+DAAbO8tLQ0Vq5cyfbt24mJiaFGjRr8+9//ZsKECTg6OpZAbUtWWloay5cvZ9euXURHR2Nvb8/jjz/OqFGj8myUCjm7x3/xxRfs27ePW7duUbt2bQYMGMCLL76IRlPxPvLJycksW7aMPXv2EBMTQ/Xq1enWrRsTJkygWrVqecpK2xUsNDSUUaNGERgYyPz58/PkSdvltWDBApYtW2Yxz8nJiZMnT5peX79+nUWLFnHkyBESExNp2LAhw4YNY9CgQbaqbqlz4MABvvnmG8LDw1GpVHh6ejJq1Ch69eqVp1xx33eyT5IVzZ8/n5UrV9K4cWOeeuopYmJi2LFjB5UrV2bDhg15Nrqs6JYtW8aCBQtwd3c3C5J0Oh1jxowhLCyMzp078/jjj3Pq1CnCwsJo2bIlq1evRqvVllDNbS81NZWhQ4dy4cIFfH19adOmDSkpKezatYuUlBRee+01xo0bB0BSUhJDhgzh8uXLdO/enfr163Pw4EHOnj1Ljx49WLRoUQm/G9tKTU1lyJAhXLx4kfbt2+Pr68vly5fZt28f7u7ubNiwAXd3d0Da7n5SU1N55plniI6Opk+fPnmCJGk7c2PHjuXAgQNMmDABlUqVJ8/e3p7x48cDOQHS4MGDSUxMpHfv3tSoUYOQkBCuXr3Kiy++yMyZM0ui+iVq9erVfPTRR1SvXp2ePXtiNBrZuXMnCQkJzJgxg9GjRwM2uu8UYRVnz55VvLy8lCFDhihZWVmm9N27dyteXl7KuHHjSrB2pUdmZqYya9YsxcvLS/Hy8lK6dOliVmbVqlWKl5eX8t///jdP+ocffqh4eXkp33zzja2qWyosWLBA8fLyUmbPnq0YjUZT+t9//6106tRJ8fHxUSIjIxVFudNGQUFBpnJ6vV6ZNGmS4uXlpezcudPm9S9J8+bNU7y8vJTFixfnSV+zZo3i5eWlvPXWW6Y0abuCzZw50/S5nT59ep48aTtznTt3VgICAu5b7pVXXlG8vLyU/fv3m9IyMjKUQYMGKd7e3kp4eHhxVrPUuXjxouLr66v07t1buXXrlin91q1bSseOHRVfX18lJSVFURTb3HcyJ8lKgoKCAJg0aVKeXo6AgADatm3L/v37iY2NLanqlQp79+6lV69erF+/nieffDLfct9//z2VKlViwoQJedKnTZuGk5MTwcHBxV3VUmX79u2oVCqmT5+e5y9Sd3d3hgwZgsFg4Ndff0Wn07FhwwZq1arF4MGDTeXs7OyYMWMGQIVru+vXr1OjRg3GjBmTJ71v374ApscSSdsVbO/evWzatImuXbua5UnbmUtISCAuLg4fH58Cy8XExLB37178/PzyfCc6ODgwffp0FEVh/fr1xV3dUmXNmjVkZ2fzwQcfUL16dVN69erVmTZtGv379yc+Pt5m950ESVZy7NgxNBoNrVu3Nsvr0KEDiqIQGhpaAjUrPX788UfS0tJ47733WL58ucUyN2/e5OrVqzRr1gxnZ+c8eU5OTjRv3pzIyEj+/vtvW1S5VHjhhRd49dVXqVy5sllebkCelpbG2bNnycjIoF27dqjVeT/adevWpX79+hw9ehSDwWCTepcGixYt4uDBg2bz2C5dugTAI488AiBtV4CEhARmzZpF69ateeGFF8zype3MnTt3DgBvb+8Cy504cQJFUejQoYNZnp+fH5UqVapwvxv79u2jZs2aFn9LBw4cyJw5c2jQoIHN7jsJkqzAYDAQGRmJh4eHxbky9evXB+Dy5cu2rlqpMnLkSPbs2cPQoUPNxuhz5f545Td/qyK25bBhw0zzF+6mKAq7d+8Gcr6MH6TtdDod169fL77KlnJJSUns3LmTadOmodFoTL2V0nb5e//990lPT+eTTz4x+zECaTtLcoOktLQ0xo0bR4cOHWjZsiXDhw/nt99+M5XLbbvc77W7aTQaatWqxfXr19HpdLapeAnL7YF77LHHiIuL45133qFTp040a9aMgQMHEhISYiprq/tOgiQrSE1NRVEUqlSpYjHf1dUVgJSUFFtWq9Rp164dLi4uBZbJbSM3NzeL+bltmZycbN3KlUHff/89p0+fpl69enTp0kXa7j7tZ4D+AAARcklEQVTWrVtH27ZtmTJlCrGxsfz3v/81/QUvbWfZL7/8ws6dO3n99dct/pCDtJ0luUHSqlWrAOjfvz9dunTh1KlTvPzyy6bpGbltkt9vR+XKlTEajaSmptqg1iUvLi4OyPlN7d+/P0eOHKFnz5707NmTS5cuMXHiRNasWQPY7r6reOsyi0F6ejpAviuuctOzsrJsVqeyKi0tDZC2vJ9t27bx0UcfodFomDt3Lvb29tJ291GtWjVefvllbt68ye7du3njjTeIi4vjxRdflLazIDY2lg8//JB27doxdOjQfMtJ25nTaDTUqVOHOXPm0LlzZ1N6REQEw4YN4+OPP6ZLly4P/NtRUXqScu+l06dP0759e5YuXYqTkxOQ03M0aNAg5s2bR9euXW1230lPkhVUqlQJgOzsbIv5uTd47v/ZIn+5bZnfl0Ju+r3zlSqS77//nunTpwMwb94809j9g96HFbXtevToweuvv868efPYunUr1apVY+7cuUREREjbWfD222+j1+v5+OOP8x0eB7nvLJk3bx579+7NEyAB+Pr6MnLkSPR6Pdu2bZPfjnvY2dmZ/j1r1qw879vT05MRI0aQnZ3Njh07bHbfSZBkBS4uLqjV6nyH03LTc7v/RP5yu07v15b3G7Yrj4xGI3PnzuWDDz7A3t6ezz//nMDAQFN+bpd9ft3LFbnt7lW7dm3Gjh0LQEhIiLTdPdatW8fvv//OjBkzqFu3boFlpe0eTtOmTQGIioq6b9slJyejUqkqTNvl/kY6OTnh6elplv/EE08AcO3aNZvddzLcZgVarZb69esTHR1NdnY29vb2efKvXbsGQOPGjUuiemVKo0aNgDttdq+K2pY6nY7p06eza9cu3Nzc+OKLL8xWf+R+qRTUdk5OTtSuXbvY61sa6HQ60+oWf39/s/zcCZ/x8fF06dIFkLbLtW3bNgBmz57N7NmzzfK3bNnCli1b6Nevn2lXaGm7HDqdjvPnz2M0GmnRooVZfkZGBpCzzL+gz6xerycmJoZHH33U4oT58qhevXpoNBr0ej2Kopj1YOb2Gjk6Otrs+06CJCtp06YNkZGRnDhxgnbt2uXJO3z4MCqVCj8/vxKqXdnh7u5OgwYNOHPmDOnp6Xm6W9PT0zl9+jQNGjSgRo0aJVhL2zIajUydOpW9e/dSt25dVq5caQom7+br64uzszNhYWEYjcY8X6zXr1/n2rVrdOzYMU+XdnlmMBgYO3YsDg4OHD582GzuQkREBAANGjSQtrtHv379aNu2rVl6dHQ0mzdvxsvLi+7du+Pj4yNtd4+0tDQGDRqEs7MzoaGhZn80Hz16FMjpUWrTpg0qlYojR44wadKkPOWOHz9OVlYWrVq1slndS5pWq6VFixYcO3aMo0ePmv2WnjlzBoAmTZrY7L6rGOGpDQwYMADIeV5PZmamKT0kJISwsDC6du2Kh4dHSVWvTBk4cCAZGRksXLgwT/qCBQtIT08vcBJpebR8+XL27t1L7dq1+f777y0GSJAzNyQwMJDr16+zevVqU7rBYGDevHlAznYCFYWjoyNPP/00qampZvfS2bNn+frrr3FyciIwMFDa7h79+/dn8uTJZv/169cPyNlyYvLkyQQEBEjb3aNq1ap07NiR1NRUlixZkifv0KFD/Pjjj3h4eNCzZ088PDzo1KkTYWFheZa3Z2Zm8tlnnwFUuO+73Pc7d+7cPNMuzp8/z7p163Bzc7PpfSfPbrOiOXPmEBQURMOGDenWrRuxsbFs376dqlWrEhwcLM9uu4e3t3e+z24bPHgwERERtG3blhYtWpie3da6dWtWrVpVYZ7dlpSUxFNPPUV6ejrdunXLdwff1q1b06FDBxISEhg4cCDR0dH861//onHjxhw6dIiIiAh69erFggULCpyEW97ExsYyZMgQoqOj8fPzo0WLFty4cYM9e/YA8Nlnn9G9e3cAabsHcOTIEV544QWzZ7dJ2+UVGRnJ0KFDiY+Pp3Xr1jRr1oyrV6+yb98+HBwc+Prrr00jC1euXGHw4MGkpKTQq1cv3N3d2bNnD5GRkYwZM4Y333yzhN+N7b311lts2rQJd3d3unfvTlpaGtu3b0ev17Nw4ULTQ71tcd9JkGRFiqIQFBTE+vXriYyMxM3NjXbt2jF16lQJkCzIL0gCTH+F7dixg/j4eDw8POjduzcvv/xyhZnECPD777+bPVLDkvHjxzNt2jQgZ6+Rzz//nP3795OSkkLdunXp378/L7zwQoUJLu+WkJDAl19+yZ49e4iLi6Ny5cq0bduW8ePHmwWd0nYFyy9IAmm7e8XGxvLFF1/w66+/cuvWLdzc3OjYsSMTJ06kYcOGecpGRkaycOFCDh8+TFZWFg0bNmTYsGEMHDiwQgWXuRRFYdOmTaxbt46//voLe3t7WrRowSuvvGI2baW47zsJkoQQQgghLJA5SUIIIYQQFkiQJIQQQghhgQRJQgghhBAWSJAkhBBCCGGBBElCCCGEEBZIkCSEEEIIYYEESUIIIYQQFkiQJIQQQghhgQRJQgghhBAWSJAkhBBCCGGBBElCCCGEEBZIkCSEEEIIYYEESUKIQtu0aRPe3t5m/z3xxBO0b9+e4cOHExwcjMFgMDt28eLFeHt7s2HDhkJdOzU1lW+++aaob0EIIfKlKekKCCHKviZNmhAQEGB6nZGRQXx8PIcOHeK9995j69atrFy5EkdHR1OZtm3bMmnSJB5//PFCXbNHjx7Y2dkxevToItdfCCEskSBJCFFkPj4+TJ482Sw9LS2Nt956i507dzJjxgwWLVpkymvXrh3t2rUr9DVv3bqFu7t7oY8XQoj7keE2IUSxcXZ2Zt68eTRs2JCdO3dy6tSpkq6SEEI8MAmShBDFytHRkVGjRgGwZcsWU7qlOUmZmZl8+umn9OnThxYtWtC6dWtGjBjBtm3bTGVy50EBxMbG4u3tzcyZM035V69eZfbs2Tz99NM0a9aM5s2b07t3bxYuXEhmZmaeunl7ezNhwgQuXLjA+PHjad26NS1atGD48OEcPnzY7L0oisIPP/zAoEGD8PPzo3379owcOZKDBw+alb106RLTp0+nU6dOPPHEE3Tr1o158+aRlJRUuIYEvv76a7y9vRk4cGCeeV56vZ7nnnsOb29vNm/eXOjzCyHykiBJCFHs2rZtC0BYWFiB5V599VVWrFhBzZo1GTZsGP/+97/5888/mTZtGuvXrwdyhvYmTZoE5PRUTZo0yTQf6vz58/Tv35+ff/6ZZs2aMXLkSAIDA4mPj2fp0qW8+eabZteMjIxk8ODBxMfH89xzz+Hv78+xY8f4//buNqSp9o8D+Het3ZVaKYoVmBnVUevFTJuRYoZETzjMMEtC7AkF61WEuVEmhTTLsNQKDJKiBzAqTbNCUytIey5N1oJKp0ZUyiof0i2v/4v+7e78nbu71z968/28cr9d53eO59V313WdbfPmzXjx4sWw69u1axfev3+P2NhYLF++HEajEZs2bcKVK1fs4xobG5GQkIDr169Do9Fg/fr1CAgIwIkTJ5CYmIju7m6X7uOGDRug0WjQ3NyMU6dO2etFRUVoampCbGws4uPjXepNRA4IIiIXXbhwQUiSJHbs2OF0XF9fn5AkSYSHh9trBQUFQpIkUVpaKoQQwmQyCUmSxPbt22XHms1mMWfOHLFkyRJZXZIkERUVJaulpaUJSZLEvXv3ZPUPHz6I0NBQERQUJD5//izrIUmSMBgMsvGFhYVCkiSRnZ1tr1VUVAhJkkRKSoqsR1tbmwgNDRURERHCarWKL1++iMjISBESEiJaWlpkfc+ePSskSRLbtm1zer+caW9vF3PnzhVqtVq0t7eLhw8fiuDgYBETEyO7LiL6dZxJIqLfTqVSAfj22P5IhBAAgFevXslmWqZOnYqrV6+ivLz8H8+TnJyM3NxcaDQaWd3b2xuzZs3C0NAQLBbLsONSU1Nlr2NiYgAAZrPZXisrKwMAZGZmwsPDw1739/eHTqfDxo0b0dfXh9raWrx//x5JSUnDntxLSkpCQEAArl275vReOOPn5wedTof+/n5kZ2cjMzMTCoUCeXl5susiol/Hp9uI6Lfr7e0F8G15bCSBgYHQaDS4f/8+oqOjERYWhoiICCxcuBBBQUE/dZ7IyEgAgMVigclkgtlshtlsRktLC549ewYAGBoakh3j6ekJLy8vWW38+PEAgMHBQXvNaDTCzc3Nvh/qRwkJCfa/m5ubAQCvX79GYWHhsLFKpRI2mw0mkwlhYWE/9X/9r9WrV6O2tha1tbUAvi0Dzp0716VeRDQyhiQi+u06OjoAfJsVcub48eMoKSlBZWUlGhoa0NDQgIMHDyIgIAB6vR7R0dFOj3/37h0MBgOuX78Om80GAJg0aRJCQ0MxadIkdHR02GesvhszZsywPgqFYljNYrHAy8vL4Xs/+vTpEwDIQowjv7KBGwCWLVtm7+9q2CIi5xiSiOi3+75he968eU7HjRs3Dunp6UhPT8fbt2/R2NiImpoaVFdXY8uWLbh27Rr8/PwcHiuEQGpqKoxGo33T98yZMzFx4kQAQGJioj2sucLd3R29vb0QQgwLSgMDAxg9ejSUSqV9tiw/Px8rVqxw+XzOdHd3Izc3F+7u7hgcHIRer0d5ebnTmToi+ve4J4mIfqvBwUGcO3cOABAXFzfiuCdPnsBgMNi/S2ny5MlYuXIlioqKsGrVKlitVjx+/HjE400mE4xGIzQaDbKyshAWFmYPSFarFa2trQAwbCbpZwUGBqKvr2/YE28AkJeXB7VajQcPHiA4OBgA0NTU5LDP0aNHcezYMYd7o35WVlYWurq6kJGRgbS0NLS3t8NgMLjcj4gcY0giot9mYGAAOp0ObW1t0Gq1Tn+CpKenByUlJSgsLJTtGxJCoLOzE4B8uU6lUtmX1IC/l80+fvwoq3/9+hX79u2zL2/9+N6/8f3R+gMHDqC/v99e7+joQHl5OTw8PKBWq7F48WJ4enrizJkzePr0qaxHVVUVDh8+jLKyMkyYMMGl67h48SKqq6sRHh6ONWvWIC0tDbNmzUJpaSnq6+td6klEjnG5jYh+mdFolG1SHhgYwNu3b3Hnzh10dXVh/vz52LNnj9MeERERWLRoEerr66HVahEREQGlUonGxkYYjUYsXboUISEh9vFTpkyB2WyGXq9HeHg44uLiEBoaikePHiEhIQELFiyA1WrF7du30draCm9vb3R1dbk8gxMfH48bN26gpqYGcXFxiIqKgtVqRVVVFXp7e1FcXAyVSgWVSoX9+/dj69atSEpKQkxMDKZNm4aXL1/i5s2bcHNzg8FgwKhRf39G/X7vUlJSnIanN2/eICcnB2PHjkVOTg4UCgX++usv5OTkYO3atdi5cycqKiqGbUQnItcwJBHRL3v+/DmeP39uf61SqeDp6YnZs2cjNjYWWq0WSqXSaY9Ro0bh0KFDOH36NCoqKnDp0iXYbDZMnz4dOp0O69atk43PysrC3r17cfnyZXR2dmLlypU4cuQICgoKcOvWLZw+fRo+Pj6YMWMG9Ho9LBYLMjIyUFdX9497oxxRKBQoKCjAmTNncPHiRZw/fx5KpRJqtRrp6en2L8wEgOjoaJw/fx7FxcW4e/cu6uvr4evrC61Wi7S0NMyYMUPWu6ioCMC3IDZSSBJCIDMzEz09PdixYwf8/f3t76nVaiQnJ+PkyZPYvXu37DfyiMh1CuHqAj0REf1f2Gw2qNVq1NXVwdfX909fDhH9F/ckERH9YZWVlfDx8YGPj8+fvhQi+gFDEhHRH2Sz2VBRUYH8/HzZPiUi+vO43EZERETkAD+2EBERETnAkERERETkAEMSERERkQMMSUREREQOMCQREREROcCQREREROQAQxIRERGRAwxJRERERA4wJBERERE5wJBERERE5ABDEhEREZEDDElEREREDjAkERERETnAkERERETkwH8AyQOSslFBmgcAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 600x450 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(4,3), dpi=150)\n",
    "plt.plot(x, sol, 'k', label='Analytical')\n",
    "plt.scatter(\n",
    "    lab.x[::10],\n",
    "    lab.species['O2'].concentration[:, -1][::10],\n",
    "    marker='x',\n",
    "    label='Numerical')\n",
    "plt.xlim([x[0], x[-1]])\n",
    "ax = plt.gca()\n",
    "ax.ticklabel_format(useOffset=False)\n",
    "ax.grid(linestyle='-', linewidth=0.2)\n",
    "plt.legend()\n",
    "plt.tight_layout()\n",
    "plt.xlim(0,60)\n",
    "plt.ylim(0,None)\n",
    "plt.xlabel('Distance, x')\n",
    "plt.ylabel('Concentration, C')\n",
    "plt.tight_layout()\n",
    "plt.savefig('analytical_transport.pdf')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}
